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Chapter 1 

INTRODUCTION 



In this work we are concerned with the one mode boson creation a' and anni- 
hilation a operators satisfying the commutation relation [a, a^] = 1 . We are 
interested in the combinatorial structures arising in the problem of normal 
ordering of a wide class of boson expressions. It shall provide us with the 
effective tools for systematic treatment of these problems. 
By normal ordering of an operator expressed through the boson a and 
operators we mean moving all annihilation operators to the right of all cre- 
ation operators with the use of commutation relation. This procedure yields 
an operator which is equivalent (in the operator sense) to the original one 
but has a different functional representation. It is of both mathematical and 
physical interest. From the mathematical point of view it allows one to rep- 
resent the operator by a function of two, in a sense 'commuting', variables. It 
is also connected with the physicist's perspective of the problem. Commonly 
used, the so called coherent state representation (see Appendix implic- 
itly requires the knowledge of the normally ordered form of the operators in 
question. Put the other way, for the normally ordered operator its coherent 
state matrix elements may be immediately read off. This representation is 
widely used e.g. in quantum optics. We mention also that calculation of 
the vacuum expectation values is much easier for operators in the normally 
ordered form. For other applications see |KS85j . 

A standard approach to the normal ordering problem is through the Wick 
theorem. It directly links the problem to combinatorics, i.e. searching for 
all possible contractions in the boson expression and then summing up the 
resulting terms. This may be efficiently used for solving problems with finite 
number of boson operators (especially when one uses the computer algebra 
packages). Although this looks very simple in that form it is not very con- 
structive in more sophisticated cases. The main disadvantage is that it does 
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not give much help in solving problems concerning operators defined through 
infinite series expansions. To do this we would have to know the underlying 
structure of the numbers involved. This still requires a lot of careful analysis 
(not accessible to computers). In this work we approach these problems us- 
ing methods of advanced combinatorial analysis |Com74j . It proves to be an 
efficient way for obtaining compact formulas for normally ordered expansion 
coefficients and then analyzing theirs properties. 

A great body of work was already put in the field. In his seminal paper 
Kat74j Jacob Katriel pointed out that the numbers which come up in the 
normal ordering problem for (a^a)" are the Stirling numbers of the second 
kind. Later on, the connection between the exponential generating function 
of the Bell polynomials and the coherent state matrix elements of e^'^^'^ was 
provided [KatOOj |Kat02j . In Chapter El we give a modern review of these 
results with special emphasis on the Dobihski relations. We also make use 
of a specific realization of the commutation relation [a, o)] = 1 in terms of 
the multiplication X and derivative D operators. It may be thought of as an 
introduction to the methods used later on in this text. Chapter E] is written 
in that spirit. We use this methodology to investigate and obtain compact 
formulas for the coefficients arising in normal ordering of a boson monomial 
(string of boson creation and annihilation operators). Then we proceed to 
the normal ordering of powers of a boson string and more generally homo- 
geneous boson polynomial {i.e. the combinations of the of the boson strings 
with the same excess of creation over annihilation operators). The numbers 
appearing in the solution generalize Stirling and Bell numbers of the second 
kind. For this problem we also supply the exponential generating functions 
which are connected with the exponentials of the operators in question. In 
each case we provide the coherent state matrix elements of the boson expres- 
sions. 

Recalling the current state of the knowledge in the field we should also men- 
tion the approach based on Lie-group methodology jWil67j . It proves useful 
for normal ordering problems for the exponentials of expressions which are 
quadratic in boson operators |Meh77j [AMTTj . 

In a series of papers by various authors |Wit75j |Mik83j |Mik85j |Kat83j some 
effort in extending these results to operators having the specific form {a} a + 
{a)y)'^ was made. In Chapter El we systematically extend this class to op- 
erators of the form {q{a))a + v{a)))'^ and e.'^i<ii°''')<^+'"i'>-'')) where q{x) and v{x) 
are arbitrary functions. This is done by the use of umbral calculus methods 
|Rom84j in finding representations of the monomiality principle {i.e. repre- 
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sentations of the Heisenberg-Weyl algebra in the space of polynomials) and 
application of the coherent state methodology. Moreover, we establish a one- 
to-one connection between this class of normal ordering problems and the 
family of Sheffer-type polynomials. 

These two classes of problems extend the current state of knowledge on the 
subject in two different directions in a quite systematic way. We believe this 
to be a considerable push forward in the normal ordering problem. We also 
emphasize the use of combinatorial methods which we give explicitly and 
which prove to be very efficient in this kind of analysis. 
We illustrate this approach by examples taken both from combinatorics and 
physics {e.g. for the hamiltonian of a generalized Kerr medium). Moreover 
we also comment on applications and extensions of the formalism. We give 
just a few of them in Chapter IHl First we observe how to extend this formal- 
ism to the case of deformed commutation relations. In general it can be done 
at the expense of introducing operator- valued Stirling numbers. Next, we 
use the fact that the use of the Dobihski relations allow us to find the solu- 
tions to the Stieltjes moment problem for the numbers arising in the normal 
ordering procedure. We use these solutions to construct and analyze new 
families of generalized coherent states |KS85j . We end by commenting on 



the specific form of the formulas obtained in Chapter El which may be used 
to derive the explicit action of generalized shift operators, called the substi- 
tution theorem. These three examples already show that the methods we 
use in solving normal ordering problems lead to many diverse applications. 
At the end this text we append complete list of publications which indicate 
some other developments on the subject. 



Chapter 2 



PRELIMINARIES 



Abstract 



In this chapter we give basic notions exploited later on in the text. We start by fix- 
ing some conventions. Next we recall the occupation number representation along 
with the creation and annihilation a operator formalism. This serves to define 
and comment on the normal ordering problem for boson operators. We end by 
pointing out a particular representation of the above in terms of the multiplication 
and derivative operators. 



2.1 Conventions 

In the mathematical literature there is always a certain freedom in making 
basic definitions. Sometimes it is confusing, though. For that reason it is 
reasonable to establish explicitly some conventions in the beginning. 
In the following by an indeterminate we primarily mean a formal variable in 
the context of formal power series (see Section^. It may be thought as a 
real or complex number whenever the analytic properties are assured. 
We frequently make use of summation and the product operations. We give 
the following conventions concerning their limits 

N N 




and 




for N <N< 



n=No 



n=No 
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for nonnegative integer k and indeterminate x. 

By the ceiling function [x] , for x real, we mean the nearest integer greater 
or equal to x. 

2.2 Occupation number representation: Boson operators 

We consider a pair of one mode boson annihilation a and creation a"!" operators 
satisfying the commutation relation 

[a,a^] = l. (2.1) 

Together with the identity operator the generators {a, a^,l} constitute the 
Heisenberg- Weyl algebra. 

The occupation number representation arises from the interpretation of a 
and as operators annihilating and creating a particle (object) in a system. 
From this point of view the Hilbert space Ti of states (sometimes called Fock 
space) is generated by the number states \n), where n = 0, 1, 2, ... count the 
number of particles (for bosons up to infinity). We assume here the existence 
of a unique vacuum state |0) such that 

a|0) = 0. (2.2) 

Then the number states {|?t-)}^o ^"6 taken as an orthonormal basis in 
H, i.e. 

{n\k) = 5n,k (2.3) 

and 

oo 

Y^\n){n\ = l. (2.4) 

n=0 

The last operator equality is called a resolution of unity which is equivalent 
to the completeness property. 

It may be deduced from Eqs. ()2.1|) and ()2.2|) that operators a and act on 
the number states as 

a \n) = y/n \n — 1), 
a^n) = v^rT+T |n + 1). 



(2.5) 
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(Indetermined phase factors may be incorporated into the states.) 
Then all states may be created from the vacuum through 

In) = -^{a^m. (2.6) 



It also follows that the number operator N counting the number of particles 
in a system and defined by 

N\n) = n\n) (2.7) 

may be represented as = a^a and satisfies the the following commutation 
relations 



[a, AT] 

[alN] = -a 



t An = _.t (2-8) 



This construction may be easily extended to the multi-boson case. 
The canonical form of the commutator Eq. (j2.ip originates from the study of 
a quantum particle in the harmonic oscillator potential and quantization of 
the electromagnetic field. Any standard textbook on Quantum Mechanics 
may serve to survey of these topics. 

We note that the commutation relations of Eqs. ()2.1|l and ()2.8j) may be easily 
extended to the deformed case |Sol94j (see also Section inHJ. 



2.3 Normal ordering 

The boson creation and annihilation a operators considered in previous 
section do not commute. This is the reason for some ambiguities in the 
definitions of the operator functions in Quantum Mechanics. To solve this 
problem one has to additionally define the order of the operators involved. 
These difficulties with operator ordering led to the definition of the normally 
ordered form of the boson operator in which all the creation operators 
stand to the left of the annihilation operators a. 

There are two well defined procedures on the boson expressions yielding a 
normally ordered form. Namely, the normal ordering M and the double dot 
: : operations. 

By normal ordering of a general function -F(a, a^) we mean M [-F(a, a^)] 
which is obtained by moving all the annihilation operators a to the right 
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using the commutation relation of Eq. (j2.1j) . We stress the fact that after this 
normal ordering procedure, the operator remains the same 
F{a, a^). It is only its functional representation which changes. 
On the contrary the double dot operation operation : F{a, a') : means the 
same procedure but without taking into account the commutation relation of 
Eq. (|2.ip . i.e. moving all annihilation operators a to the right as if they 
commute with the creation operators a^. We emphasize that in general 
F{a, a^) 7^ : F{a, a^) : . The equality holds only for operators which are 
already in normal form {e.g. Af [F{a, a^)'j = : Af [F{ 



a, a 



Using these two operations we say that the normal ordering problem for 
F{a, a^) is solved if we are able to find an operator G{a, a"!") for which the 
following equality is satisfied 

F{a,a^)=Ar[F{a,a^)]=:G{a,a^): . (2.9) 

This normally ordered form is especially useful in the coherent state repre- 
sentation widely used in quantum optics (see Appendix Also calculation 
of the vacuum expectation values in quantum field theory is immediate when- 
ever this form is known. 

Here is an example of the above ordering procedures 

aa^aaa^a — > (a^)^a^ + 4 a^a^ + 2 

\a,an=l ^ ^ ■' 



- to the left a - to the right 
aa^aaa^a '■ — ' > a'^a'^aaa 

a, at - commute 
(like numbers) 

Another simple illustration is the ordering of the product a'^'(a^)' which is in 
the so called anti-normal form [i.e. all annihilation operators stand to the 
left of creation operators). The double dot operation readily gives 

: a\a^y : = (at)'a^ 

while the normal ordering procedure Af requires some exercise in the use of 
Eq. ()2.1|) yielding (proof by induction) 



' = AT [a^a^Y] = J2 (^) ^- {ci^y~''a''-P. (2.10) 
p=o VP/ 
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These examples explicitly show that these two procedures furnish completely 
different results (except for the operators which are already in normal form). 
There is also a 'practical' difference in their use. That is while the application 
of the double dot operation : : is almost immediate, for the normal ordering 
procedure H certain skill in commuting operators a and is needed. 
A standard approach to the problem is by the Wick theorem. It reduces the 
normal ordering procedure M to the double dot operation on the sum over 
all possible contractions {contraction means removal of a pair of annihilation 
and creation operators in the expression such that a precedes a^). Here is an 
example 

aa^ aaa'^ aaa'^ — : { all contractions } : 

= : aa^aaa^aaa^ : 

+ : aa^aaa^a ^ + aa^aaa^ ^a + aa}a (^a^aa ^ + 
aa} <^aa}aa f£^+ ^a^aaa^aa ifi + aa}a ^ ^^aaa^ + 
aa^ ^a ^'^aaa'^+ ^a'^aa ^'^aaa'^+ ^ ^'^aaa^ aaa'^ : 

+ : aa^ a ^, ^'^a ^ ^ + aa!^ ^a i^^a ^ ^a^aa (fia (fi ^ 
f£ ^'^aaa^a ^ + aa^a f£a + aa^ ^a f£a + 
fia^aa ifi ^a + ^ fi^aaa^ f£a + aa^ ^ ^ ^^aa + 
^a^a (^aa <fi \ ^ fi^a ^a^aa + aa^ ^ ^'^aa + 
(^a ^aa ^ (^aci'aa (^a^a ^ ^aa ^ + 
(^a) ^,a, aa ifi ^ (fi a ^ f£^aaa^+ ^ (f} (fiaaa} : 

+ : (pa ^ (pa ^ p + ^ p ^a pa ^ p + ^ pa ^ p fla p + 
^ P fia p ^a p+ ^ p fi fi paa p+ fi p ^ ^ paa p : 

= (a^) V + 9 (a^)^a^ + 18 a^a^ + 6 al 

One can easily see that the number of contractions may be quite big. This 
difficulty for polynomial expressions may be overcome by using modern com- 
puter algebra systems. Nevertheless, for nontrivial functions (having infinite 
expansions) the problem remains open. Also it does not provide the ana- 
lytic formulas for the coefficients of the normally ordered terms in the final 
expression. A systematic treatment of a large class of such problems is the 
subject of this work. 
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At the end of this Section we recall some formulas connected with operator 
reordering. The first one is the exponential mapping formula, sometimes 
called the Hausdorff transform, which for any well defined function -F(a, a^) 
yields 

e^'^Ffa, at)e-^" = F(a,a^ + x), 

t . (2-11) 

e-^''V(a, a"f)e^"^ = F{a + x,a^). 

It can be used to derive the following commutators 
[a,F(a,at)] = ^F(a,at), 

(2.12) 

[at,F(a,at)] = -d_F{a,a)). 

The proofs may be found in any book on Quantum Mechanics, e.g. fLou64j . 
Also a well known property of the Heisenberg-Weyl algebra of Eq. (j2.1|) is a 
disentangling formula 

gA(a+at) ^ e^'/^gAat^Aa ^ ^-X^ /2 ^Xa ^Xa^ ^ (2.13) 

which may serve as an example of the normal ordering procedure 



e 



A(a+at) ^ 



gA(a+at) 



= e^'/2 : e^('^+'^') : . (2.14) 



This type of expressions exploits the Lie structure of the algebra and uses a 
simplified form of Baker- Campbell-Hausdorff formula. For this and other dis- 
entangling properties of the exponential operators, see |Wil67j |Wit75j |MSI98j |Das96j . 
Finally, we must mention other ordering procedures also used in physics, like 
the anti-normal or the Weyl (symmetric) ordered form. We note that there 
exist translation formulas between these expressions, see e.g. |CG69j(5T79j . 



2.4 X and D representation 

The choice of the representation of the algebra may be used to simplify the 
calculations. In the following we apply this, although, we note that with 
some effort one could also manage using solely the operator properties of a 
and . 

There are some common choices of the Heisenberg-Weyl algebra representa- 
tion in Quantum Mechanics. One may take e.g. a pair of hermitian operators 
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X and P = i-^ acting on the dense subspace of square integrable functions 
or the infinite matrix representation of Eq. ()2.5p defined in the Fock space. 
This choice is connected with a particular interpretation intimately connected 
with the quantum mechanical problem to be solved. 

Here we choose the simplest possible representation of the commutation re- 
lation of Eq. (|2.H) which acts in the space of (formal) polynomials. We make 
the identification 

" n (2-15) 

where X and D are formal multiplication and derivative (D = ^) operators, 
respectively. They are defined by their action on monomials 



n+l 



Xx"' = X , 



(2.16) 



(It also defines the action of these operators on the formal power series.) 
Note that the commutation relation of Eq. (j2.H) remains the same, i.e. 



[D,X] = 1. (2.17) 
For more details see Appendix lEl and Chapter El 

This choice of representation is most appropriate, as we note that the prob- 
lem we shall be concerned with has a purely algebraic background. We are 
interested in the reordering of operators, and that only depends on the al- 
gebraic properties of the commutator of Eqs. ()2.1|) or ()2.17p . We emphasize 
that the conjugacy property of the operators or the scalar product do not 
play any role in that problem. To get rid of these unnecessary constructions 
we choose the representation of the commutator of Eq. ()2.1|) in the space of 
polynomials defined by Eq. ()2.16|) where these properties do not play primary 
role, however possible to implement (see the Bargmann-Segal representation 
|Bar61j|Seg63| ). We benefit from this choice by the resulting increased sim- 
plicity of the calculations. 



Chapter 3 



GENERIC EXAMPLE. STIRLING AND 
BELL NUMBERS. 

Abstract 



We define Stirling and Bell numbers as solutions to the normal ordering problem 
for powers and the exponential of the number operator N = a}a , where a and 
are the annihilation and creation operators. All their combinatorial properties are 
derived using only this definition. Coherent state matrix elements of A^" and e'^^ 
are shown to be the Bell polynomials and their exponential generating function, 
respectively. 



This Chapter may be thought of as a simple introduction to some methods 
and problems encountered later on in the text. We give a simple example 
on which all the essential techniques of the proceeding sections are used. 
Although Stirling and Bell numbers have a well established purely combina- 
torial origin |(]om74j |Rlo84j |Wil94j [r;KP94j we choose here a more physical 
approach. We define and investigate Stirling and Bell numbers as solutions 
to the normal ordering problem KatTi] [KatOOj |Kat02j . 

Consider the number operator N = a^a and search for the normally 
ordered form of its n-th power (iteration). It can be written as 

n 

{a^af = J2S{n,k){a^)''a'' (3.1) 

k=l 

where the integers S{n, k) are the so called Stirling numbers. 
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One may also define the so called Bell polynomials 

n 

B{n,x) = ^S{n,k)x^ (3.2) 



k=l 

and Bell numbers 

n 

B{n) = B{n,l) = J2S{n,k) (3.3) 

k=l 

For convenience we apply the conventions 

S{n, 0) = 6nfl and S{n, k) = ior k > n (3.4) 

and 

B{0) = B{0,x) = 1. (3.5) 

(It makes the calculations easier, e.g. one may ignore the summation limits 
when changing their order.) 

In the following we are concerned with the properties of these Stirling and 
Bell numbers. 

First we state the recurrence relation for Stirling numbers 

S{n + l,k) = kS{n,k) + S{n,k-l), (3.6) 

with initial conditions as in Eq. (j3.4|) . The proof 9by induction) can be de- 
duced from the equalities 

71+1 

S{n + 1, k) (a"^)^a^ = (a"^a)"^^ = a'^a (a'^a)'' 

k=l 

= J2 Sin, k) a^aia^fa"" ^ S{n, k) {a^f{a^a + k)a^ 

k=l k=l 

n+1 n 

= J2 S{n, k-1) {a^)''a'' + ^ S{n, k) {a^fa"" 



k=2 k=l 
n+1 



^ J2 iS{n, k-l) + k S{n, k)) {a^fa^. 



k=l 
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We shall not make make use of this recurrence relation, preferring the sim- 
pler analysis based on the Dobihski relation. To this end we proceed to the 
essential step which we shall extensively use later on, i.e. change of represen- 
tation. Eq. ()3.1|) rewritten in the X and D representation (see Section 12. 4p 
takes the form 

n 

{XDf = J2 S{n, k) X^DK (3.7) 

k=l 

With this trick we shall obtain all the properties of Stirling and Bell numbers. 
We first act with this equation on the monomial x™. This gives for m integer 



Y,S{n,k)mK 

k=l 

where x- = x ■ {x — 1) ■ ... ■ {x — k + 1) is the falling factorial {x- = 1). 
Observing that a (non zero) polynomial can have only a finite set of zeros 
justifies the generalization 



n 



X 

k=l 



J2S{n,k)xK (3.8) 



This equation can be interpreted as a change of basis in the space of poly- 
nomials. This gives an interpretation of Stirling numbers as the connection 
coefficients between two bases {a;"}^o {^~}^o- 
Now we act with Eq. (j3.7j) on the exponential function e^. We obtain 



oo I, 
.X 



fc=o ■ fc=l 
Recalling the definition of the Bell polynomials Eq. (j3.2j) we get 

oo ^ 

5(n,x)=e-^^r|-, (3.9) 

k=0 

which is a celebrated Dobinski formula |Wil94j . It is usually stated for Bell 
numbers in the form 

k=l 
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Note that both series are convergent by the d'Alembert criterion. 
The most striking property of the Dobihski formula is the fact that integer 
numbers B{n) or polynomials B{n, x) can be represented as nontrivial infinite 
sums. Here we only mention that this remarkable property provides also 
solutions to the moment problem (see Section 16. 2j) 

In the following we shall exploit the Dobinski formula Eq. (j3.9|) to investigate 
further properties of Stirling and Bell numbers. 

Applying the Cauchy multiplication rule Eq. ()B.3|) to Eq. ()3.9|) and comparing 
coefficients we obtain explicit expression for S{n, k) 



Next, we define the exponential generating function of the polynomials B{n, x) 
(see Appendix in} as 

G{X,x) = Y,B{n,x) — , (3.12) 

n=0 

which contains all the information about the Bell polynomials. Substituting 
Eq. ()3.9|) into Eq. ()3.12|) . changing the summation order and then identifying 
the expansions of the exponential functions we obtain 



OOOO h \ ^ OOi.OO -.^ 



n=0 k=0 k=0 n=0 



oo ^ 

^ k\ 



k=0 

Finally we may write the exponential generating function G{X,x) in the 
compact form 

G(A,a;) =e"(^'-i). (3.13) 

Sometimes in applications the following exponential generating function of 
the Stirling numbers S{n, k) is used 

T,SMi^=^^^- (3-14) 

n=k 
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It can be derived by comparing expansions in x of Eqs. (j3.12j) and (j3.13|) 
G{\,x) = 1 + ^ K]] A;): 



oo / n 
n=l \k=l 

k=l \n=k 



G(A,x) = e^(^'-i) = 1 + ^ ^^-^x^ 

fc=i 

Bell polynomials share an interesting property called the Sheffer identity 
(note resemblance to the binomial identity) 

Bin, x + y) = J2C'^ ^(^' - ^' ^) (3-15) 

fc=0 ^ ^ 

It is the consequence of the following equahties 
n! 



(x+y)(e^-l) ^ gi^(e^-l)gy(e^-l) 



n=0 



E5(n,a; + ?/)— = e 



A" -^^^ .A" 



n\ — ' n\ 

n=0 n=0 

n=0 \fc=0 ^ ^ / ■ 

By this identity and also by the characteristic exponential generating function 
of Eq. ()3.13|) Bell polynomials B{n,x) are found to be of Sheffer-type (see 
Appendix IB]). 

Differentiating the exponential generating function of Eqs. (j3.12j) and (|3.13|) 
we have (see Eq. flB.5|) ) 



V Bin + l,x)^ = ^G(A, x) = eV(^'-i) = V B(n, x)^ 
^-^ n\ OA ^-^ n\ 

n=0 n=0 



which by the Cauchy product rule Eq. ()B.3jl yields the recurrence relation for 
the Bell polynomials 

B{n\\,x) = ^{^\B{k,x). (3.16) 
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And consequently for the Bell numbers we have B{n + 1) = Y12=o iT)^^^^' 
By the same token applied to the exponential generating function of Stirling 
numbers S{n, k) of Eq. ()3.14|) the recurrence relation of Eq. ()3.6p can be also 
derived. 

Finally, using any of the derived properties of Stirling or Bell numbers one 
can easily calculate them explicitly. Here are some of them 



S{n,k), l<k<n B{n) 



n = 


1 


1 














1 


n = 


2 


1 


1 












2 


n = 


3 


1 


3 


1 










5 


n = 


4 


1 


7 


6 


1 








15 


n = 


5 


1 


15 


25 


10 


1 






52 


n = 


6 


1 


31 


90 


65 


15 


1 




203 


n = 


7 


1 


63 


301 


350 


140 


21 


1 


877 


n = 


8 


1 


127 


966 


1701 


1050 


266 


28 1 


4140 



Now we come back to normal ordering. Using the properties of coherent 
states (see Appendix we conclude from Eqs. ()3.1|) and ()3.2|) that diagonal 
coherent state matrix elements generate Bell polynomials |KatOOj 

{z\{a^a)''\z) = B{n,\z\^). (3.17) 

Moreover, expanding the exponential e^"^^" and taking the diagonal coherent 
state matrix element we have 

{z\e''^'^\z) = J2{^\{a^<^r\z)- = Y,B{n, \z\')-. 

n=0 n=0 

We see that the diagonal coherent state matrix elements of e'*"*^" yield the 
exponential generating function of the Bell polynomials (see Eqs. (j3.12j) and 

(EIISl)) 

(^le^'^'^l^) = el"l'(^'-^\ (3.18) 
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This result allows us to read off the normally ordered form of e^"" " (see 
Appendix |XJ 



These are more or less all properties of Stirling and Bell numbers. Note 
that they were defined here as solutions to the normal ordering problem for 
powers of the number operator = a^a. The above analysis relied firmly on 
this definition. 

On the other hand these numbers are well known in combinatorial analysis 
|Com74j |Rio84j |Wil94j [GKP94j where are called Stirling numbers of the sec- 
ond kind and are usually explored using only recurrence relation Eq. (j3.6p . 
Here is their original interpretation is in terms of partitions of the set. 

• Stirling numbers S{n, k) count the number of ways of putting n different 
objects into k identical containers (none left empty). 

• Bell numbers B{n) count the number of ways of putting n different 
objects into n identical containers (some may be left empty). 




X e 



(3.19) 



Some other pictorial representations can be also given e.g. in terms of graphs 
[MBP05^ or rook numbers \N iw7'A\ [SPB+Oi] | Var()4j . 



Chapter 4 



NORMAL ORDERING OF BOSON 
EXPRESSIONS 



Abstract 



We solve the normal ordering problem for expressions in boson creation a' and 
annihilation a operators in the form of a string (a^)^'"a**f . . . {a^Y'^a^'^ (a^)^'^a*^. We 
are especially concerned with its powers (iterations). Next we extend the solution 
to iterated homogeneous boson polynomials, i.e. powers of the operators which are 
sums of boson strings of the same excess. The numbers obtained in the solutions 
are generalizations of Stirling and Bell numbers. Recurrence relations, closed- 
form expressions (including Dobihski-type relations) and generating formulas are 
derived. Normal ordering of the exponentials of the the aforementioned operators 
are also treated. Some special cases including Kerr-type hamiltonians are analyzed 
in detail. 



4.1 Introduction 

In this chapter we consider expressions in boson creation and annihilation 
a operators (see Section [221) • search for the normally ordered form and 
show effective ways of finding combinatorial numbers arising in that problem. 
These numbers generalize Stirling and Bell numbers (see Chapter O]). 
The first class of expressions treated in this Chapter are so called boson 
strings or boson monomials, see Section 14.21 Here is an example 

aa'' aaa^ a' a'' aa'' aaaaaaa'' aa^ aaa'' aa^ , 
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which we compactly denote as 

with some nonnegative integers and Sm- By an excess of a string we call 
the difference between the number of creation and annihilation operators, 
=i{^m ~ Sm)- In the following we assume that the boson strings 
are of nonnegative excess (for negative excess see Section H^H)) • 
Boson strings (monomials) can be extended further to polynomials. We treat 
homogeneous boson polynomials which are combinations of boson monomials 
of the same excess (possibly with some coefficients). Here is an example 

aaa) aa) a) aaa) a) a) + a)aa)a)a + a\ 

which using methods of Section 14.21 can be written in the form 

N 
k=No 

with appropriate c?, A^, Nq and a^'s. 

Next, in Sections 14.31 and 14.41 we give recipes on how to approach the prob- 
lem of normal ordering of iterated boson strings and homogeneous boson 
polynomials, i. e. their n-th powers, like 

{a) aa^ aaa"^ a'^ a'^ aa^ aaaaaaa"^ aa'^ aaa^ a'^ aa^^ 

and 

These considerations serve to calculate in Section 031 the exponential gener- 
ating functions of generalized Stirling and Bell numbers. This in turn gives 
the solution to the normal ordering problem of the exponential of the afore- 
mentioned boson operators. 

In Section we comment on operators with negative excess. 

Finally in the last Section 14.71 we work out some examples which include 

generalized Kerr-type hamiltonian. 

We refer to |MBPn5j (BPS+?)5j |SDB+n4j |Varn4j for interpretation of consid- 
ered combinatorial structures in terms of graphs and rook polynomials. 
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4.2 General boson strings 

In this section we define the generahzation of ordinary Bell and Stirling 
numbers which arise in the solution of the general normal ordering problem 
for a boson string (monomial) |MBP05j [WitOSj . Given two sequences of 
nonnegative integers r = (ri, r2, . . . , tm) and s = (si, S2, ■ ■ ■ , sm) we define 
the operator 

^r,s = {a^y^'a'^' . . . {a^y^a'^a^Y'a'K (4.1) 

We let Sr,s{k) be the nonnegative integers appearing in the normally ordered 
expansion 

S1+S2H hsji/ 

H,,s = {a^r'' ^r,s(fc) (at)'a^ (4.2) 

k=si 

where dn = X]m=i(^m ~ "^m); ''^ = 1, M. Here we assume that the (overall) 
excess of a string is nonnegative > (for negative excess see Section 
Observe that the r.h.s. of Eq. (j4.2j) is already normally ordered and consti- 
tutes a homogeneous boson polynomial being the solution to the normal or- 
dering problem. We shall provide explicit formulas for the coefficients Sr^s{k)- 
This is a step further than the Wick's theorem which is nonconstructive in 
that respect. 

We call S'r,s(^) generalized Stirling number. The generalized Bell polynomial 
is defined as 

\-sm 

5r,s(x)= Yl S^Mx", (4.3) 

k=si 

and the generalized Bell number is the sum 

S1+S2H hSM 

5r,s = 5r,s(l) = Yl '^r,s(A;). (4.4) 

k=si 

Note that sequence of numbers Sr,s{k) is defined here for si < k < si + 
S2 + ■ ■ ■ + Sm- Initial terms may vanish (depending on the structure of r 
and s), all the next are positive integers and the last one is equal to one, 

5'r,s(Sl + S2 H h Sm) = 1- 

For convenience we apply the convention 

Sr,s{k) =0 for /c < si or > si + S2 + ■ ■ ■ + sm- (4.5) 
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We introduce the notation r l+l tm+i = (^i, r2, . . . , tm, '^m+i) and s l±) sm+i = 
(si, S2, . . . , sm, sm+i) and state the recurrence relation satisfied by general- 
ized Stirhng numbers Sr,s{k) 



5r«.M^-.,s«.M.-.(^) = E ''''''' (^^'^ + ^ - ^^'-^^ ^r,s(A: - j)- (4.6) 
i=o V J / 

One can deduce the derivation of Eg . (14.61) from the following equalities 

siH hSA/+i 

k=s\ 

— -'-'rl±)rM+i,sl+)SM+l — / " -'-'r,s 



5^ (at)'^a' 

fc=Sl 

fc=si 

SlH hSA/ SAf+l / s 



fc=Sl 

SlH VsM 

k=s\ 

SlH hSM SAf+l 



fc=si j=0 
SA/+1 ^ V SlH f-SM+l-j 

j,-=0 ^ ^ ' fe=Si+SM+l-j 
S1---+SA,/ SAf+i , , 

k=si j=0 V / 

I hSA/ SM + 1 y s 



fc=si j=0 
SlH hSA/ sm+1 



[a 

k=si j=0 
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The problem stated in Eq. (j4.2j) can be also formulated in terms of multipli- 
cation X and derivative D operators (see Section |22 



X''W' ...X''WX''W = X'^" J2 Sr,s{k)X^D^ (4.7) 

k=s\ 

Acting with the r.h.s. of Eq. (j4.7j) on the monomial we get x'^'^" 'Yli'k=s^^ Sr,s{k)l- 
On the other hand action of the l.h.s. on x' gives nm=i('^ni-i + I'^^x^^'^^' . 
Equating these two results (see Eq. ()4.7p ) we have 

M s-l+S2^ hSAf 

m=l k=si 

By invoking the fact that the only polynomial with an infinite number of 
zeros is the zero polynomial we justify the generalization 

M S1+S2A \-sm 

ll{d^.i + x)^= Yl Sr,s{k)xK (4.8) 



m=l 



k=si 



This allows us to interpret the numbers Sr^s{k) as the expansion coefficients 
of the polynomial {nm=i('^"^-i + x)—} in the basis {x-}^q. 
Now, replacing monomials in the above considerations with the exponential 
we conclude that 



S1+S2H hSM 



M 



m=l 



X 



k=si k=si 

With these two observations together and Eq. ()4.7|) we arrive at the Dobinski- 
type relation for the generalized Bell polynomial 

M 



BrJx) 



k=si 



,m=l 



X 
Id 



(4.9) 



Observe that the d'Alembert criterion assures the convergence of the series. 
Direct multiplication of series in Eq. ()4.9|) using the Cauchy rule of Eq. ()B.3|) 
gives the explicit formula for generalized Stirling numbers Sr,s{k) 



1 /k\ ^ 

7=0 ^"'^ m=l 



(4.10) 
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We finally return to normal ordering and observe that the coherent state 
matrix element (see Appendix^ of the boson string yields the generalized 
Bell polynomial (see Eqs. ()4.1|) and dOl)) 

{z\H,Az) = (z*)^«5r,s(kr). (4.11) 
4.3 Iterated homogeneous polynomials in boson operators 

Consider an operator of the form 

TV 
k=No 

where the a^'s are constant coefficients (for convenience we take a^o and 
Oat 7^ 0) and the nonnegative integer d is the excess of the polynomial. Note 
that this is a homogeneous polynomial in boson operators. It is already in a 
normal form. If the starting homogeneous boson polynomial is not normally 
ordered then the methods of Section provide proper tools to put it in the 
shape of Eq. ()4.12j) . For boson polynomials of negative degree see Section 1^?^ 
Suppose we want to calculate the normally ordered form of (if^)". It can be 
written as 

nN 
k=No 

with Sci{n,k) to be determined. These coefficients generalize Stirling num- 
bers (see Eq. fl3.1|) ). In the same manner as in Eqs. ()3.2|) and ()3.3p we define 

generalized Bell polynomials 

nN 

Bi{n,x)= Y,St{n^)x\ (4.14) 

k=No 

and generalized Bell numbers 

nN 

Biin) = Bi{n,l) = Si{n,k). (4.15) 

k=No 
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Note that S^{n, k), n, x) and B^{n) depend on the set of a^s defined in 
Eq. ()4.12j) . For convenience we make the conventions 

si{o,o) = 1, 

Si{n,k) = for k > nN, (4.16) 
s£{n, fc) = for A; < A^'o and n > 0, 

and 

51(0) = 51(0, x) = l. (4.17) 

In the following we show how to calculate generalized Stirling and Bell num- 
bers explicitly. First we state the recurrence relation 

N I . . 

Si{n + l,k) = J2(^i J2( ){nd + k-l+p)^Si{n,k-l+p), (4.18) 

/=Aro p=0 

with initial conditions as in Eq. ()4.16p . 

It can be deduced from the following equalities 



{n+l)N 



fat 



1113. t 



){n+i)d J2 Si{n + 1, k) (at)'=a^ ^ (^1)"+' = H^{H, 

k=No 

N nN 

)^ ^az(at)'a' {a^'' J] ^^n, A;) (at)'=a'= 

l=No k=No 
nN N 

J2 Si{n, k) J2 [a^y^'ci' (at)"''+^a*^ 

1 l=No 

nN N / ' /7\ 

Siin, k) ^ a, {ar^^ E ( J + 



k=No l=No 



k=No l=No \p=0 

nN N ' //\ 

fc=Aro l=No P=0 

i=Aro p=0 k=No+l~p 

■ , J n. nN N ' //\ 

(^t)(n+i)d ^ ^ ^ / M (^rf + ^ _ / + p)P 5^^^^ ^ _ ^ ^ p)(at)'=a^ 

k=Nol=No p=0 
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Eq. ()4.13p rewritten in the X and D representation takes the form 



nN 



k=No 



(4.19) 



where 



N 



H^{D,X)=X'^ ^afcX'^D^ 

k=No 



(4.20) 



We shall act with both sides of the Eq. (j4.19|) on the monomials xK The r.h.s. 
yields 



nN nN 
A:=Aro k=No 



■nd 



On the other hand iterating the action of Eq. (j4.2(Jj) on monomials we have 

n N 



.1=1 k=No 



X 



l+nd 



Then when substituted in Eq. fl4.19p we get 

n N nN 



n E 

i=l k=No 



m-= Y.^'c.{n,k)A 

k=No 



Recalling the fact that a nonzero polynomial can have only a finite set of 
zeros we obtain 



n N nN 

J] ^ ttfc (x + (z - l)d)^ = J2 '^"("' ^) 

i=l k=No k=No 



(4.21) 



This provides interpretation of the generalized Stirling numbers S^{n,k) 
as the connection coefficients between two sets of polynomials {x-}^q and 
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Now acting with both sides of Eq. ()4.19|) on the exponential function we get 

nN 
oo 



1=0 



' n N 

n 5^ ttfc + - 1)^^)- 

.1=1 fc=Afo 



7[ 



Taking these equations together we arrive at the Dobinski-type relation 



n N 



1=0 \_i=l k=No 



X 
1. 



(4.22) 



Note that by the d'Alembert criterion this series is convergent. 

Now multiplying the series in Eq. ()4.22p and using the Cauchy multiplication 

rule of Eq. ()B.3|) we get the explicit expression for generalized Stirling numbers 



1 ^ /k\ " ^ 

j=0 ^ i=l l=No 



ai (j + (z - l)dy-. 



(4.23) 



We end this section observing that the diagonal coherent state matrix ele- 
ments (see Appendix of (Ha)"' generate generalized Bell polynomials (see 
Eq.dl^D) 

{z\{H'S\z) = {zT'Bi{n,\zf). (4.24) 

Other properties of generalized Stirling and Bell numbers will be investigated 
in subsequent Sections. 



4.4 Iterated boson string 

In Section 1^21 it was indicated that every boson string can be put in normally 
ordered form which is a homogeneous polynomial in a and . (Note, the 
reverse statement is evidently not true, i.e. 'most' homogeneous polynomials 
in a and can not be written as a boson string.) This means that the results 
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of Section mav be applied to calculate the normal form of iterated boson 
string. We take an operator 

^r,s = {a^Y"a'^' . . . {a^Y^a'^ {a^Y'a'' (4.25) 

where r = (ri, r2, . . . , r^) and s = (si, S2, ■ ■ ■ , sm) are fixed integer vectors 
and search for its n-th power normal form 

n{si+S2A I-sm) 

(ffr,s)" = {a^Y' E ^r,s(n, k){a^Ya', (4.26) 

k=si 

where d = Y.f=i{'^m - s^)- 

Then the procedure is straightforward. In the first step we use results of 
Section lOl to obtain 

S1+S2H hSM 

Hr,. = {a^Y E S,,,{kYa^Ya\ 

k=si 

One can use any of the formulas in the Eqs. ()4.6p . ()4.9|) or ()4.10j) . 

Once the numbers Sr,s{k) are calculated they play the role of the coefficients 

in the homogeneous polynomial in and a which was treated in preceding 

Section. Taking = Sj.,s{k), Nq = si and = Si + S2 + ... + s^/ in Eq. (j4.12|) 

we arrive at the normal form using methods of Section 

We note that analogous scheme can be applied to a homogeneous boson 

polynomial which is not originally in the normally ordered form. 

4.5 Generating functions 

In the preceding Sections we have considered generalized Stirling numbers as 
solutions to the normal ordering problem. We have given recurrence relations 
and closed form expressions including Dobihski-type formulas. Another con- 
venient way to describe them is through their generating functions. We define 
exponential generating functions of generalized Bell polynomials B'^{n,x) as 

°° A" 
GiiX,x) = J2Biin,x)- (4.27) 

n=0 

They are usually formal power series. Because B'^{n) grows too rapidly with 
n these series are divergent (except in the case = for k > 1 treated in 
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Chapters El and EI) • In Section UTI we give a detailed discussion and show 
how to improve the convergence by the use of hypergeometric generating 
functions. 

In spite of this 'inconvenience', knowing that formal power series can be 
rigorously handled (see Appendix^, in the following we give some useful 
(formal) expressions. 

Substituting the Dobihski-type formula of Eq. ()4.22j) into definition Eq. ()4.27p 



oo oo 



n=0 1=0 

which changing the summation order yields 



n N 

n E + - 1)^)- 

i=l k=No 



oo ; OO 

X 



' n N 

n 5^ ttfc (/ + (z - i)d)^ 

1=1 k=No 



(4.28) 



1=0 n=0 

Using the Cauchy multiplication rule of Eq. ()B.3|) we get the expansion in x: 



Gi{X,x) 



oo J 
X-' 



3=0 •' 1=0 



n=0 



' n N 

n E + - 1)^)- 

i=l k=No 



(4.29) 



which by comparison with 

oo / n 

G^(A,x) = l + $^ 5^5^(n,A;) 



X 



n=l \k=l 



A" 



oo / oo 



1+ E 

k=No \n=k 



A'^ 



X 



gives the exponential generating function of the generalized Stirling numbers 
S^{n, k) in the form (compare with Eq. ()4.23|) ) 



A 



n=k 



1=0 



n N 



n E + - 1)^)- 



n=0 \_i=l k=No 



A^ 



(4.30) 



We note that for specific cases Eqs.dOHl, and (|Om) simphfy a lot, 

see Section 14.71 for some examples. 
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This kind of generating functions are connected with the diagonal coherent 
state elements of the operator e^^^^ . Using observation Eq. ()4.24j) we get 

n=0 ■ n=0 

which yields 

{z\e^''^z) = Gi{{z*Y\\z\^) . (4.32) 

This in turn lets us write (see Appendix ^ the normally ordered expression 

^XHi ^ . G<'^{{a^yX,a^a) : . (4.33) 

Observe that the exponential generating functions (^^(Ajx) in general are 
not analytical around A = (they converge only in the case = for 
k > 1, see Chapters El and El) • Nevertheless, e.g. for negative values of XaN 
and d = the expressions in Eqs. ()4.28|) or ()4.29|) converge and may be used 
for explicit calculations. The convergence properties can be handled with 
the d'Alembert or Cauchy criterion. 

We also note that the number state matrix elements of e^^« are finite (the 
operator is well defined on the dense subset generated by the number states). 
It can be seen and explicitly calculated from Eqs. ()4.33|) . ()4.28|) or ()4.29|) . 

4.6 Negative excess 

We have considered so far the generalized Stirling and Bell numbers aris- 
ing in the normal ordering problem of boson expressions and their powers 
(iterations) with nonnegative excess. Actually when the excess is negative, 
i.e. there are more annihilation a than creation a''" operators in a string, the 
problem is dual and does not lead to new integer sequences. 
To see this we first take a string defined by two vectors r = (ri, r2, . . . , r^) 
and s = (si, S2, ■ ■ ■ , sm) with a negative (overall) excess dM = X]m=i(^m ~ 
Sm) < 0. The normal ordering procedure extends the definition of the num- 
bers Sr,s{k) for the case of negative excess through 

ri+r2A hTM 

k=rM 

(4.34) 
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Definition of Br^s{x) and Bj- g is analogous to Eqs. (j4.3|) and (j4.4j) : 

ri+r2^ hrM 

5r,s(x)= Yl •^r,s(fc)x^ (4.35) 

k=rM 

and 

ri+r2A hVM 

k=rM 

Note that the limits in the above equations are different from those in Section 
By taking the hermitian conjugate we have 

ri+r2^ hrM , i n j i , 

k=rM 

+ JJTl / ri+r2+-+rM 



k=rM 



T\+r2-\ l-r-M 



k=rM 



where s = {sm, ■ ■ ■ ,S2, si) and r = {tm, • • • , ^"2, ''"i) correspond to the string 
ifs,r with positive overall excess equal to —dM > 0. 
This leads to the symmetry property 

SrAk) = 5s,f(fc) (4.37) 

as well as 

Br,s{n^x) = Bs^{n,x) and Br^s{n) = Bs^{n) (4.38) 
for any integer vectors r and s. 

Formulas in Section may be used to calculate the generalized Stirling and 
Bell numbers in this case if the symmetry properties Eqs. ()4.37j) and ()4.38j] 
are taken into account. We conclude by writing explicitly the coherent state 
matrix element of a string for du < 0: 

{z\{a^Y''a''' . . . {a^Y'a''{a^Y"a''\z) = B,4\z\^)z~'^^ . (4.39) 
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Now we proceed to homogeneous polynomials of Section 14.31 In the case of 
the negative excess we define 

N 
k=No 

with some constant coefficients a^'s (for convenience we take a^o 7^ and 
«iv 7^ 0) and nonnegative integer d. The excess of this homogeneous boson 
polynomial is —d. We search for the numbers S~'^{n, k) defined by the normal 
form of the n-th power of Eq. ()4.40|) as 

nN 

{H-JT = E ^"'(^, k) {a^fa' a^" . (4.41) 

fc=Afo 

Generalized Bell polynomials B~'^{n, x) and Bell numbers B~'^{n) are defined 
in the usual manner. 

In the same way, taking the hermitian conjugate, we have 

^ S~An,k) {a^a' a^'^ ™ ((i^.^)")')' = {{{H^yTy 



k=No 
nN 

k=No 

(Note that (5^.)* = Si , see any formula of Section Oil 
Consequently the following symmetry property holds true 

S-\n,k) = St{n,k) (4.42) 

and consequently 

B-\n,x) = Bi{n,x) and B-\n) = B^^n) (4.43) 
for any vector ex and integer d. 

Again any of the formulas of Section 14.31 may be used in calculations for the 
case of negative excess. Coherent state matrix element of Eq. ()4.4()j) takes the 
form 

{z\{H-JY\z) = B-J{n,\z\')z'. {AAA) 
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4.7 Special cases 

The purpose of this section is to illustrate the above formalism with some ex- 
amples. First we consider iteration of a simple string of the form (a"'")'" a* and 
work out some special cases in detail. Next we proceed to expressions involv- 
ing homogeneous boson polynomials with excess zero d = 0. This provides 
the solution to the normal ordering of a generalized Kerr-type hamiltonian. 

4.7.1 Case {{a^Ya^y 

Here we consider the boson string in the form 

Hr,s = {a^Ya' (4.45) 

for r > s (the symmetry properties provide the opposite case, see Section 

EM- 

It corresponds in the previous notation to Hr^s with (r, s) = (r, s) or with 
d = r — s and = 1 (zero otherwise). 

The n-th power in the normally ordered form defines generalized Stirling 
numbers Sr,s{n,k) as 

ns 

{HtS = (at)"(^-«)^5,,,(n,fc)(at)'=a^ (4.46) 

k=s 

Consequently we define the generalized Bell polynomials Bj.^s{n,x) and gen- 
eralized Bell numbers Br^s{n) as 

ns 

Br,s{n, x)=J2 Sr,s{n, k)x^ (4-47) 

k=s 

and 

ns 

Br,s{n) = Br,s{n, 1) = 5^ Sr,s{n, k). (4.48) 

k=s 

The following convention is assumed: 
5.,s(0,0) = 1, 

Sr,s{n,k) = for > ns, (4.49) 
Sr^s{n, k) = for k < s and n > 0, 
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and 



(4.50) 



As pointed out these numbers and polynomials are special cases of those 
considered in Section lO} i.e. 



for d = r — s and = 1 (zero otherwise). 

In the following the formalism of Section 14.31 is applied to investigate the 
numbers in this case. By this we demonstrate that the formulas simplify 
considerably when the special case is considered. The discussion also raises 
some new points concerning the connection with special functions. 
Below we choose the way of increasing simplicity, i.e. we start with r > s 
and then restrict our attention to r = s and s = 1. For r = 1 and s = 1 we 
end up with conventional Stirling and Bell numbers (see Chapter IH}. Where 
possible we just state the results without comment. 

Before proceeding to the program as sketched, we state the essential results 
which are the same for each case. We define the exponential generating 
function as 



Sr,s{n,k) = S^^\n,k) 



oo 





ra=0 



The coherent state matrix elements can be read off as 




(4.52) 



and 




(4.53) 



This provides the normally ordered forms 




(4.54) 
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and 

e^^" = : a,s(A(at)^-^ata) : . (4.55) 

We shall see that usually the generating function of Eq. ()4.51|1 is formal and 
Eq. ()4.5Hj) is not analytical around A = 0. No matter of the convergence sub- 
tleties operator equations Eqs. ()4.54|l and ()4.55|1 hold true in the occupation 
number representation. 



r>s 

Recurrence relation 



Sr,s{n + l,k) = J2 (^) (^(^ -s) + k-r + p)^ 5,,,(n, k - s + p) (4.56) 



with initial conditions as in Eq. ()4.49p . 

Connection property 



Hix + (j - l)(r - s)r = Yl SrAn, k)xK 
i=i 



(4.57) 



k=s 



The Dobihski-type relation 



k=s j=l 

-s)'("~^)e^^^ 

k=0 



II ^ r—s 

11 r(i + ^) 

.1=1 \ r—s' 



X 
Id 



(4.58) 
(4.59) 



Explicit expression 



p=s ^-^^ 7 = 1 



ip+ij-l){r-s)f. (4.60) 



The "non-diagonal" generalized Bell numbers Br^s{n) can always be expressed 
as special values of generalized hypergeometric functions pFg , e.g. the se- 
ries i?2r,r(^) can be written down in a compact form using the confluent 
hypergeometric function of Kummer: 



(rn)! 



e ■ r\ 



iFi(ra + l,r + 1; 1) 



(4.61) 
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and still more general family of sequences has the form (p, r = 1, 2 . . 



{p{n-l)+])\ 

(pjV- 



n 

rFr{pn + 1, pn + 1 + p, . . . , pn + 1 + p{r - 1); 1 + p, 1 + 2p, . . . , 1 + rp; 1) 



etc... 

The exponential generating function takes the form 



oo I CO 
1=0 ■ n=0 



n(^ + (^-l)(r-.))^ 



i=l 



n=l 



n(/ + (^-l)(r-.))^ 



A" 



For s > 1 it is purely formal series (not analytical around A = 0). 
A well-defined and convergent procedure for such sequences is to consider 
what we call hypergeometric generating functions, which are the exponen- 
tial generating function for the ratios Br^s/{n\y, where t is an appropriately 
chosen integer. A case in point is the series i?3 2(^) which may be written 
explicitly from Eq. ()4.58|) as: 



oo 



{n + k)\{n + k + l)\ 



e— A;!(A; + l)!(fc + 2)! 



Its hypergeometric generating function 6*3,2 (A) is then: 



G3,2(A) = 



n=0 



A" 



Similarly for 6*4,2 (A) one obtains: 



, - + 1; 1;4A 

2 ' 2 
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Eq. (j4.61|) implies more generally: 

'B2r,r{ny 



G'2r,r(A) — ^ 



n=0 



iTJ Efclo (fcTT)!2^i I + ^' ^' ^'^)' 

j_ ■spoo 1 p (k+1 k±2 fc+3. 1 i.97\N 
2!e /^fc=0 (fe+2)!3-^2 3 , 3 , 3 , J-, J-, ^ ' , 



r = 2, 
r = 3, 



3bEr=o7IT3)T4^3(Y,...,^;l,l,l;256A), r = 4... 



etc... 



See jJPSOlj for other instances where this type of hypergeometric generating 
functions appear. 



r=s 

Recurrence relation 



Srrin + 1, = ](k - r + p)^ SrJn, k - r + p) 



with initial conditions as in Eq. ()4.49|) . 
Connection property 

nr 

k=r 

The Dobihski-type relation has the form: 

Br,r{n,x) = e-'^^k^— 

'{k + r)\ 



k=s 
00 r / 

E 

k=0 



k\ 



IT' 



Explicit expression 



p=r ^ ' 



(4.62) 



(4.63) 



(4.64) 
(4.65) 



(4.66) 
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Using Eq. ()4.66p we can find the following exponential generating function of 

E ^) = h-^y - 1) ' (4-67) 

n=\k/r] p=r 

We refer to Section 14.7.21 for considerations of the exponential generating 
functions. 

s=l 

Recurrence relation 

Sr,iin + l,k) = {n{r - 1) + k - r + k) + Sr,i{n, k - 1) (4.68) 

with initial conditions as in Eq. ()4.49|) . 
Connection property 

n ns 

Y[{x + (j - l)(r - s)) = J2 Sr,i{n, k)xK (4.69) 

i=l k=s 

The Dobihski-type relation has the form: 

oo n f, 

Br^n, x) = e-^ 5] n + - - • (4-70) 

A:=l j=l 

Explicit expression 

The generalized Bell numbers Br^i{n) can always be expressed as a combina- 
tion of r — 1 different hypergeometric functions of type iF^._i(. . . ; x), each of 
them evaluated at the same value of argument x = {r — 1)^^^ ; here are some 
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lowest order cases: 



B2,i{n) 
B3,i{n) 

54,i(n) 



-iFi(n + l;2;l) = (n-l)!4^2i(-l), 

e 



(4.72) 



)n— 1 



2r(n + i; 



iF2(n+l;l,l;\] + n^F^ ( n + 1; ^, 2; ^ 



2' 2' 2' 4 



2' '4 



-ii^3ln+-,-,-,-,-) + 



2e 



TT 



3r(n + 



2 2 4 5 1 \ 



iFAn + ] + nhFs U + 1; -, -, 2; — 



3 3' 3' 3 27 J 



4 5 



3 3 '27 



etc... 



In Eq. fj4.72p Lm\y) is the associated Laguerre polynomial. 

In this case the exponential generating function of Eq. ()4.51|) converges (see 

also Chapter EI) 



-^l-(r-l)A 



(4.73) 



The exponential generating function of Sr,i{n, k) takes the form 



Y] --rSr,i{n,k) 



n=\k/r~\ 



(1 - (r - l)x) - 1 



(4.74) 



See also Chapter for detailed discussion of this case. 
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We end this section by writing down some triangles of generalized Stirling 
and Bell numbers, as defined by Eqs. fj4.46|) and ()4.48|) . 











5'i,i(n, /c) 


l<k<n 




n 


= 1 


1 








1 


n 


= 2 


1 


1 






2 


n 


= 3 


1 


3 


1 




5 


n 


= 4 


1 


7 


6 1 




15 


n 


= 5 


1 


15 


25 10 


1 


52 


n 


= 6 


1 


31 


90 65 


15 1 


203 









S2,l{ 


n,k), 1 < k < n 


B2,l 


n 


= 1 


1 






1 


n 


= 2 


2 


1 




3 


n 


= 3 


6 


6 


1 


13 


n 


= 4 


24 


36 


12 1 


73 


n 


= 5 


120 


240 


120 20 1 


501 


n 


= 6 


720 


1800 


1200 300 30 1 


4051 



r=3, s=l 









S3,i{n,k) 


l<k<n 


-83,1 (" 


n 


= 1 


1 






1 


n 


= 2 


3 


1 




4 


n 


= 3 


15 


9 


1 


25 


n 


= 4 


105 


87 


18 1 


211 


n 


= 5 


945 


975 


285 30 1 


2236 


n 


= 6 


10395 


12645 


4680 705 45 1 


28471 



r=2, s=2 



S'2,2(n,fc), 2 < fc < 2n 



-62,2 (n) 



1 

2 

4 

8 

16 

32 



4 1 

32 38 12 1 

208 652 576 188 24 1 

1280 9080 16944 12052 3840 580 

7744 116656 412800 540080 322848 98292 



40 

16000 



1 

1390 



60 



1 
7 

87 
1657 
43833 
1515903 



=3, s=2 



S3,2(n,k), 2 < fc < 2n 



B3,2(n) 



1 
6 
72 

1440 

43200 

1814400 



6 

168 

5760 

259200 



1 

96 

6120 
424800 



15120000 34776000 



18 

2520 

285120 

33566400 



1 

456 
92520 



36 

15600 



16304400 4379760 



1 

1380 
682200 



60 1 
62400 3270 



90 



1 

13 
355 
16333 
1121881 
106708921 



r=3, s=3 



Sa.sCn.fe), 3 < fe < 3n 



S3,3(n) 



1 

6 

36 

216 

1296 



18 9 

540 1242 

13608 94284 

330480 6148872 



1 

882 

186876 

28245672 



243 27 1 

149580 56808 11025 
49658508 41392620 18428400 



1107 54 1 
4691412 706833 63375 



3285 90 1 



1 

34 
2971 
513559 
149670844 
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4.7.2 Generalized Kerr Hamiltonian 

In this section we consider the case of the homogeneous boson polynomial 
with excess zero {d = 0). This illustrates the formalism of Sections 14.31 and 
14.51 on the example which constitutes the solution to the normal ordering 
problem for the generalized Kerr medium. 

The Kerr medium is described by the hamiltonian H = a^a^aa, |MW95j . By 
generalization we mean the hamiltonian in the form 

N 

k=No 

Observe that this is exactly the form of the operators (with excess d = 0) 
considered in Section lOl i.e. 

^ Hi (4.76) 

To simplify the notation we skip the index = in S'^{n, k), B^{n, x), B^{n) 
and C^IX, x). 

We are now ready to write down the solution to the normal ordering problem 
in terms of generalized Stirling numbers 

nN 

{H^T= Sc.in,k) (at)^a^ (4.77) 

k=No 

generalized Bell polynomials 

nN 

B^{n,x)= J2SUn,k)x', (4.78) 

k=No 

and Bell numbers 

nN 
fc=7Vo 

Recurrence relation 

TV I . . 

Sc.{n + l,k)=Y(^i J2[ jik~l+Pr-S^{n,k-l + p), (4.80) 

l=No p=0 
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with initial conditions as in Eq. (j4.16|) . 
Connection property 

N \ " nN 

\k=No 



J2 Sc.{n, k) xK 

k=No 



The Dobihski-type relation has the form: 



N 



1=0 \k=No 



Explicit expression 



N 
\l=No 



Exponential generating function 



{n,x) 



n=0 



nl 



n 



n=0 



rX 



E 

m=0 



E 

.1=0 



X 

ml 



(4.81) 



(4.82) 



(4.83) 



(4.84) 
(4.85) 
(4.86) 



Again, the exponential generating function of Eq. ()4.84|) is formal around 
A = (for > 1). Although we note that by the d'Alembert criterion for 
Aoat < the the series of Eq. (j4.85|) converges. 
The generating function of the Stirling numbers is 



-1 



(4.87) 



n=k 1=0 

Now, returning to the coherent state representation we have 



{z\e^"-\z) =G'„(A, l^n. 



(4.88) 
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This provides the normally ordered form of the exponential 

e^^« = : G„(A,a^a) : . (4.89) 

These considerations may serve as a tool in investigations in quantum optics 
whenever the coherent state representation is needed, e.g. they allow one 
to calculate the Husimi functions or other phase space pictures for thermal 
states jBHPSOBj . 



Chapter 5 



MONOMIALITY PRINCIPLE AND 
NORMAL ORDERING 

Abstract 



We solve the boson normal ordering problem for (^q{a^)a + v{a^)^^ with arbitrary 
functions q{x) and v{x). This consequently provides the solution for the exponen- 
tial gHii<^'')<^+'"iO' )) generalizing the shift operator. In the course of these considera- 
tions we define and explore the monomiality principle and find its representations. 
We exploit the properties of Sheffer-type polynomials which constitute the inher- 
ent structure of this problem. In the end we give some examples illustrating the 
utility of the method. 



5.1 Introduction 

In Chapter E] we treated the normal ordering problem of powers and exponen- 
tials of boson strings and homogeneous polynomials. Here we shall extend 
these results in a very particular direction. We consider operators linear in 
the annihilation a or creation operators. More specifically, we consider 
operators which, say for linearity in a, have the form 

g(a^)a + t>(a^) 

where q{x) and v{x) are arbitrary functions. Passage to operators linear in 

is immediate through conjugacy operation. 
We shall find the normally ordered form of the n-tli power (iteration) 



(g(a^)a + f (a^)) 
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and then the exponential 

gA(q(at)a+^;(at))_ 

This is a far reaching generahzation of the results of |Mik83j |Mik85j |Kat83j 
where a special case of the operator a) a + was considered. 
In this approach we use methods which are based on the monomiality prin- 
ciple [BDHP05J. First, using the methods of umbral calculus, we find a wide 
class of representations of the canonical commutation relation Eq. (j2.17|) in 
the space of polynomials. This establishes the link with Sheffer-type polyno- 
mials. Next the specific matrix elements of the above operators are derived 
and then, with the help of coherent state theory, extended to the general 
form. Finally we obtain the normally ordered expression for these operators. 
It turns out that the exponential generating functions in the case of linear 
dependence on the annihilation (or creation) operator are of Sheffer-type, 
and that assures their convergence. 

In the end we give some examples with special emphasis on their Sheffer-type 
origin. We also refer to Section for other application of derived formulas. 

5.2 Monomiality principle 

Here we introduce the concept of monomiality which arises from the action of 
the multiplication and derivative operators on monomials. Next we provide 
a wide class of representations of that property in the framework of Sheffer- 
type polynomials. Finally we establish the correspondence to the occupation 
number representation. 

5.2.1 Definition and general properties 

Let us consider the Heisenberg-Weyl algebra satisfying the commutation re- 
lation 

[P,M] = 1. (5.1) 

In Section 12.41 we have already mentioned that a convenient representation 
of Eq. (j5.ip is the derivative D and multiplication X representation defined 
in the space of polynomials. Action of these operators on the monomials is 
given by Eq. ()2.16p . Here we extend this framework. 
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Suppose one wants to find the representations of Eq. (j5.1j) sucli tfiat tlie action 
of M and P on certain polynomials s„(rc) is analogous to the action of X 
and D on monomials. More specifically one searches for M and P and their 
associated polynomials Sn[x) (of degree n, n = 0, 1, 2, ...) which satisfy 



The rule embodied in Eq. fl5.2|l is called the monomiality principle. The poly- 
nomials Sn{x) are then called quasi-monomials with respect to operators M 
and P. These operators can be immediately recognized as raising and low- 
ering operators acting on the s„(x)'s. 

The definition Eq. ()5.2j) implies some general properties. First the operators 
M and P obviously satisfy Eq. ()5.1|) . Further consequence of Eq. ()5.2|) is the 
eigenproperty of MP, i.e. 



The polynomials s„(x) may be obtained through the action of M" on sq 



MSn{x) = Sn+l{x), 
PSn{x) = n Sn-l{x). 



(5.2) 



MPSn{x) = nSn{x). 



(5.3) 



(5.4) 



and consequently the exponential generating function of s„(x)'s is 



oo 



A 




nl 



(5.5) 



n=0 



Also, if we write the quasimonomial Sn{x) explicitly as 



71 




(5.6) 



fe=0 



then 




(5.7) 



Several types of such polynomial sequences were studied recently using this 
monomiality principle |D()TV97j |(4DM97j |Dat99j |DS(:01j \CesOO\ . 
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5.2.2 Monomiality principle representations: Sheffer-type 
polynomials 

Here we show that if Sn{x) are of Sheffer-type then it is possible to give 
exphcit representations of M and P. Conversely, if M = M{X, D) and 
P = P{D) then Sn{x) of Eq. ()5.2|) are necessarily of Sheffer-type. 
Properties of Sheffer-type polynomials are naturally handled within the so 
called umhral calculus |Rom84j |RR78j |Huc98j (see AppendixlH]). Here we put 
special emphasis on their ladder structure. Suppose we have a polynomial 
sequence s„(x), n = 0, 1,2, ... (s„(a:) being a polynomial of degree n). It is 
called of a Sheffer A- type zero |She39j . fRai65j (which we shall call Sheffer- 
type) if there exists a function f{x) such that 

f{D)sn{x) = nsn-i{x). (5.8) 

Operator f{D) plays the role of the lowering operator. This characterization 
is not unique, i.e. there are many Sheffer-type sequences s„(x) satisfying 
Eq. ()5.8|l for given f{x). We can further classify them by postulating the ex- 
istence of the associated raising operator. A general theorem |R,om84j [CheOSj 
states that a polynomial sequence s„(a;) satisfying the monomiality principle 
Eq. ()5.2|) with an operator P given as a function of the derivative operator 
only P = P{D) is uniquely determined by two functions f{x) and g{x) such 
that /(O) = 0, / (0) 7^ and g{0) ^ 0. The exponential generating function 
of Sn{x) is then equal to 

^ g{f 1(A)) 

and their associated raising and lowering operators of Eq. ()5.2|) are given by 

P = f{D), 



M 



X 9'{D) 



1 



(5.10) 



f'(D) ■ 



Observe the important fact that X enters M only linearly. Note also that 
the order of X and D in M{X, D) matters. 

The above also holds true for f{x) and g{x) which are formal power series. 
By direct calculation one may check that any pair M, P from Eq.dEini) auto- 
matically satisfies Eq. ()5.H) . The detailed proof can be found in e.g. |R,om84j . 
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Here are some examples we have obtained of representations of the monomi- 
ahty principle Eq. ()5.2|) and their associated Sheffer-type polynomials: 

a) M{X, D) = 2X - D, P{D) = \D, 
Sn{,x) = Hn{x) - Hermite polynomials; 
G{X,x) = e^^^-^'. 

b) M(X, D) = -XD^ + (2X -1)D-X-1, P{D) = - ^'^^-^ D"^, 
Sn{x) = n\Ln{x) - where Ln{x) are Laguerre polynomials; 

c) M(X, D) = Xj^, P{D) = + D, 
Sn{.x) = Pn{x) - Bessel polynomials jGro78j : 

d) M(X, D) = X(l + D), P{D) = ln(l + D), 
Sn{x) = Bn{x) - Bell polynomials; 

G(A,x) =e^(^'-i). 

e) M(X, D) = Xe~^, P{D) = - 1, 

Sn{x) = X- - the lower factorial polynomials |TS95j : 
G{X,x) = e^'i°(i+^). 

f) M(X, D) = {X- tan(D)) cos\D), P{D) = tan(L)), 
Sn{x) = Rn{x) - Hahn polynomials |Ben91j : 

g) M(X, D) = X^^D, P{D) = Wl{D), 
where Wl{x) is the Lambert W function |CGH+96j : 
Sn{x) = In{x) - the idempotent polynomials |Com74j : 

G{X,x) =e^^^\ 
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5.2.3 Monomiality vs Fock space representations 

We have already called operators M and P satisfying Eq. ()5.2j) the rising 
and lowering operators. Indeed, their action rises and lowers index n of 
the quasimonomial s„(x). This resembles the property of creation and 
annihilation a operators in the Fock space (see Section I2.2|) given by 

til ^ ^^".V (5-11) 

These relations are almost the same as Eq. (|5.2|) . There is a difference in 
coefficients. To make them analogous it is convenient to redefine the number 
states \n) as 



\n) = Vn\ \n). (5.12) 

(Note that fo) = |0)). 

Then the creation and annihilation operators act as 

= ^J!iJ^^ (5.13) 
a^|n) = |n + 1). 

Now this exactly mirrors the the relation in Eq. ()5.2|) . So, we make the 
correspondence 

P < — > a 

M < — > (5.14) 
Sn{x) < > \n) , ra = 0,1,2,... . 

We note that this identification is purely algebraic, i.e. we are concerned 
here only with the commutation relation Eqs. ()5.1|) . ()2.17p or ()2.1|) . We do 
not impose the scalar product in the space of polynomials nor consider the 
conjugacy properties of the operators. These properties are irrelevant for our 
proceeding discussion. We note only that they may be rigorously introduced, 
see e.g. ^Rom84j. 



5.3 Normal ordering via monomiality 



In this section we shall exploit the correspondence of Section lB.2.3l to obtain 
the normally ordered expression of powers and exponential of the operators 
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a^q{a)+v{a) and (by the conjugacy property) q{a^)a+v{a^). To this end shall 
apply the results of Section 15.2.21 to calculate some specific coherent state 
matrix elements of operators in question and then through the exponential 
mapping property we will extend it to a general matrix element. In conclusion 
we shall also comment on other forms of linear dependence on a or . 
We use the correspondence of Section 15.2.31 cast in the simplest form for 
M = X, P = D and s„(x) = x", i.e. 

D < > a 



X < — > (5.15) 
< > \n) , n = 0,1,2,... . 

Then we recall the representation Eq. ()5.10|) of operators M and P in terms 
of X and D. Applying the correspondence of Eq. ()5.15|) it takes the form 



P(a) = /(a 
M(a,at) = 



at - ^ 

9{«) 



1 



(5.16) 



/'(«) ■ 



RecaUing Eas.(l01). (IKl)|l and ^J^i we get 

n 

[M(a,at)]"|0) = 5^s„,fc (at)». (5.17) 

k=0 

In the coherent state representation it yields 

{z\ [M(a, at)]" |0) = s„(z*)(z|0). (5.18) 
Exponentiating M(a, at) and using Eq. (j5.9|) we obtain 

^^|gAAf(a,at)|Q^ = ^^^^^^ 6^*^"' ^ (;2|0) . (5.19) 

By the same token one obtains closed form expressions for the following 
matrix elements ( |/) is the l-th number state, / = 0, 1, 2, ...) 

{z\ [M(a, at)]" 10 = ±^s^^i{z*){z\0), (5.20) 



and when property Eq. ()B.5|) is applied 
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Observe that in both Eqs. (j5.18p and (j5.2(jp we obtain Sheffer-type polyno- 
mials (modulus coherent states overlapping factor {z\0)). Also Eqs. ()5.19p 
and ()5.21|) reveal that property through the Sheffer-type generating func- 
tion. This connection will be explored in detail in Section 15.4.21 
The result of Eq. (j5.19|) can be further extended to the general matrix element 
^^|gAM(a,at)|y^^ To this end recall Ea. (|X3|l and write 



Next, using the exponential mapping property Eq. ()2.11|) we arrive at 



' I (;2|e r ^^(^^Wj f'(a+z') |0). 



Now we are almost ready to apply Eg. ()5.19j) to evaluate the matrix element 
on the r.h.s. of the above equation. Before doing so we have to appropriately 
redefine functions f{x) and g{x) in the following way {z' - fixed) 

fix) ^ f(x) = f{x + z')-f{z'), 
g{x) g{x) = g{x + z')/g{z'). 

Then /(O) = 0, /'(O) 7^ and ^(0) = 1 as required by Sheffer property for 
f{x) and g{x). Observe that these conditions are not fulfilled by f{x + z') 
and g{x + z'). This step imposes (analytical) constraints on z', i.e. it is valid 
whenever f'{z') 7^ (although, we note that for formal power series approach 
this does not present any difficulty). Now we can write 



'(^|0). 



^(/-HA)) 

By going back to the initial functions f{x) and g{x) this readily gives the 
final result 

/ I AM(a,at)| /\ _ gj^') „z*\,f-HX+,f(z'))-z'],\J\ 
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where {z\z') = e^*^' 2'^ 2 1^1^ is the coherent states overlapping factor (see 
Eq.dm). 

To obtain the normally ordered form of e'*'*^^"''*^) we apply the crucial property 
of the coherent state representation of Eqs. ()A.ll|) and ()A.12|) . Then Eq. ()5.22p 
provides the central result 

AAf(a,at) _ . a^rH\+f{a))~a] M . 

9{f-\X + f{a)))- ■ ^'-'^^ 

Let us point out again that appears linearly in M(a,a^"), see Eq. (l5.16|) . 
We also note that the constraints for functions /(x) and g{x), i.e. f{0) = 0, 
/'(O) 7^ and g{0) 7^ 0, play no important role. For convenience (simplicity) 
we put 

1 



g'jx) 1 

g{x) f'{x)' 



and define 



T(A,x) = /-i(A + /(x)), 
9{T{X,x)) 

This allows us to rewrite the main normal ordering formula of Eq. ()5.2Hj] as 

^X{aUia)+v{a)) ^ . ^at [T(A,a)-a] ^^^^^^ . ^^ ^A) 

where the functions T(A, x) and G{X, x) fulfill the following differential equa- 
tions 

^I^ = q{T{X,x)) , nO,x)=x, (5.25) 

dG{X,x 



v{T{X,x))-G{X,x) , G{0,x) = l. (5.26) 

In the coherent state representation it takes the form 

(^'|eA(atq{a)+.(a)) _ ^^/|^^ ^z'*[T{X,z)-z]q^^^ ^5_27) 
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We conclude by making a comment on other possible forms of linear depen- 
dence on a or a"''. 

By hermitian conjugation of Eq. ()5.24|) we obtain the expression for the nor- 
mal form of e^('''-°^^""'"^'-'^'') . This amounts to the formula 

with the same differential equations Eqs. ()5.25p and ()5.26|) for functions T(A, x) 
and G{X,x). In the coherent state representation it yields 

^^/|gA(g(at)a+.(a)) _ ^^^^ ^/*)g[T(A,.'*)-.'*]. ^^ ^9) 

We also note that all other operators linearly dependent on a or a"!" may 
be written in just considered forms with the use of Eq. (j2.12p . i.e. ag(a^) + 
f (a"'") = g(a"'")a + g'(a"'') + f (a"'') and q{a)a^ + v{a) = a"l"g(a) + q'{a) + v{a). 
Observe that from the analytical point of view certain limitations on the 
domains of z, z' and A should be put in some specific cases (locally around 
zero all the formulas hold true). Also we point out the fact that functions 
q{x) and v{x) (or equivalently f{x) and g{x)) may be taken as the formal 
power series. Then one stays on the ground of formal power expansions. 

In the end we refer to Sect ion where these results are applied to derive 
the substitution formula. We note that the reverse process, i.e. derivation 
of the normally ordered form from the substitution theorem, is also possible, 
see [BHP^ . 

5.4 Sheffer-type polynomials and normal ordering: 
Examples 

We now proceed to examples. We will put special emphasis on their Sheffer- 
type origin. 

5.4.1 Examples 

We start with enumerating some examples of the evaluation of the coherent 
state matrix elements of Eqs. ()5.18|) and ()5.27p . We choose the M{a,a^)^s as 
in the hst a) - g) in Section 15.2.21 

a) {z\{—a + 2a"l')"|0) = Hn{z*){z\0), Hermite polynomials; 

^^|gA(-a+2at)|^,^ = e^^^'*-''^-^' {z\z') . 
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b) {z\ [-ata + (2at - l)a - at + l]" |0) = n\Ln^i{z*){z\Q) , 
Laguerre polynomials; 



c) {z 

d) {z 

e) {z 

f) {z 

g) {z 
{z 



(atj^)" |0) = Pn{z*){z\Q), Bessel polynomials; 

(a+a + at)"|0) = 5„(2*)(z|0), Bell polynomials; 
(a"''e~'^)"|0) = {z*)-l^z\Qi) , the lower factorial polynomials; 



\{a) — tan(a)) cos^(a)]" |0) = Rn{z*){z\Q) ^ Halm polynomials; 

oA(at-tan(a)) cos2(a) |^/\ _ cos[arctan(A+tan(2'))] ^z* [arctanfA tanf^'D-z'] /^UA 
I ' cos{z') \ I '■ 



" WL(a) " 



|0) = In{z*){z\Q) , the idempotent polynomials; 



Aat 



l + Wj (a) , , W^ 



Note that for = we obtain the exponential generating functions of ap- 
propriate polynomials multiplied by the coherent states overlapping factor 
(^10), see Eq.(ISI22)). 

These examples show how the Sheffer-type polynomials and their exponential 
generating functions arise in the coherent state representation. This generic 
structure is the consequence of Eqs. (j5.18|) and (j5.22j) or in general Eqs. (j5.27|) 
or (j5.29p and it will be investigated in more detail now. Afterwords we shall 
provide more examples of combinatorial origin. 



5.4.2 ShefFer polynomials and normal ordering 

First recall the definition of the family of Sheffer-type polynomials s„(z) 
defined (see Section \^ through the exponential generating function as 

G(A, x) = Y, ■'^niz)— = A{X) e^^W (5.30) 

n=0 
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where functions A{\) and B{X) satisfy: A{0) ^ and 5(0) = 0, 5'(0) 0. 
Returning to normal ordering, recall that the coherent state expectation value 
of Eq. ()5.28|) is given by Eq. ()5.29|) . When one fixes z' and takes A and z as 
indeterminates, then the r.h.s. of Eq. ()5.29p may be read off as an exponential 
generating function of Sheffer-type polynomials defined by Eq. (|5.3U|) . The 
correspondence is given by 

A{\)=g{\z!*), (5.31) 
5(A) = [T(A,/*) -2'*]. (5.32) 

This allows us to make the statement that the coherent state expectation 
value {z'\...\z) of the operator exp [A(g(a"'')a + f (a^))] for any fixed z' yields 
(up to the overlapping factor (^'l^)) the exponential generating function of 
a certain sequence of Sheffer-type polynomials in the variable z given by 
Eqs. ()5.31|l and ()5.32|1 . The above construction establishes the connection be- 
tween the coherent state representation of the operator exp [A(q'(a''')a + f (a^))] 

and a family of Sheffer-type polynomials s'n'"^ (z) related to q and v through 



^^/|gA[,(at)a+.(at)]|^^ = {z'\z) 1 + V s(f'-)(^)^ , (5.33) 




where explicitly (again for z' fixed): 

s(f'")(^) = {z'\zY^{z'\ {q{a^)a^v{a^)Y\z). (5.34) 

We observe that Eq. ()5.34|) is an extension of the seminal formula of J. Katriel 
|Kat74j . [Kat00j where v{x) = and q{x) = x. The Sheffer-type polynomials 
are in this case Bell polynomials expressible through the Stirling numbers of 
the second kind (see Section |2I). 

Having established relations leading from the normal ordering problem to 
Sheffer-type polynomials we may consider the reverse approach. Indeed, it 
turns out that for any Sheffer-type sequence generated by A{X) and -B(A) one 
can find functions q{x) and v{x) such that the coherent state expectation 
value (z'lexp [X{q{a'')a + v{a^))~\ \z) results in a corresponding exponential 
generating function of Eq. (j5.3Uj) in indeterminates z and A (up to the over- 
lapping factor {z'\z) and z' fixed). Appropriate formulas can be derived from 
Eqs.dESU) and (IHIT^ by substitution into Eqs.^^ and 

q{x) = B'{B~\x - z'*)), (5.35) 
A'(B~\x - z'*)) 
^ ^ A{B-^{x - z'*)) ^ ' 
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One can check that this choice of q{x) and v{x), if inserted into Eqs. (j5.'25|) 
and ()5.26|) . results in 

r(A,x) = B{X + B~\x - z*)) + z'*, (5.37) 

^^^'"^^ - A{B-\x-z'*)) ' ^^-^^^ 

which reproduce 

^^/|gA[g(at)a+.(at)]|^^ = {z' \z) A{\)e'^'^^\ (5.39) 

The result summarized in Eqs. ()5.3H) and ()5.32|) and in their 'dual' forms 
Eqs. fj5.35|) - ()5.38p provide us with a considerable flexibility in conceiving and 
analyzing a large number of examples. 



5.4.3 Combinatorial examples 



In this section we will work out examples illustrating how the exponential 
generating function G{\) = J2'^=o^n^ of certain combinatorial sequences 
{an)'^=o appear naturally in the context of boson normal ordering. To this 
end we shall assume specific forms of q{x) and v{x) thus specifying the oper- 
ator that we exponentiate. We then give solutions to Eqs. ()5.25p and ()5.26p 
and subsequently through Eqs. ()5.31|) and ()5.32|) we write the exponential 
generating function of combinatorial sequences whose interpretation will be 
given. 



a) Choose q{x) = x"", r > 1 (integer), v{x) = (which implies g{X, x) = 1). 
Then T(A, x) = x[l — A(r — l)x^~"'^] ^""^ . This gives 



,A(a"l')'-a 



= : exp 



(1 - A(r- l)(at)^-i)' 



1 a 



as the normally ordered form. Now we take 2;' = 1 in Eqs. ()5.31|) and 
()5.32p and from Eq. ()5.39p we obtain 



exp 



(1 - A(r - 1))~ 



- 1 
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which for z = 1 generates the following sequences: 



r = 2: a„ = 1,1,3,13,73,501,4051,... 
r = 3: a„ = 1,1,4,25,211,2236,28471, 



etc. 



These sequences enumerate r-ary forests |Slo05j |Sta99j jFS05j . 

b) For q{x) = a;ln(ex) and v{x) = (implying g{X,x) = 1) we have 



T(A,x) 



■^a;^". This corresponds to 



Af 



^Aat ln(eat)a 



= : exp 



a^y - 1 a 



whose coherent state matrix element with z' = 1 is equal to 



(l|;z)-i(l|e^"""('^"')'^|;2) 



exp 



z e 



which for 2; = 1 generates an = 1, 1, 3, 12, 60, 385, 2471, ... enumerating 
partitions of partitions |Sta99j . |Slo05j . |FS05j . 

The following two examples refer to the reverse procedure, see Eqs. fj5.35|) - 
fl5.38|) . We choose first a Sheffer-type exponential generating function and 
deduce q{x) and v{x) associated with it. 

c) A{X) = B{\) = A, see Eq. fl5.3()j) . This exponential generating 

function for 2; = 1 counts the number of arrangements a„ = n\ Ylk=o h ~ 
1, 2, 5, 65, 326, 1957, ... of the set of n elements jCom74j . The solutions 
of Eqs. ()5.35p and ()5.36|) are: q{x) 
bosons it corresponds to 



1 and v(x) = In terms of 



AT 



A aH ^^ 

e V 2-at, 



,Aa 



it -A 



it - A 



d) For A{\) = 1 and B{X) = 1 — Vl — 2A one gets the exponential gen- 
erating function of the Bessel polynomials |(Tro78j . For z = 1 they 
enumerate special paths on a lattice |Pit99j . The corresponding se- 
quence is a„ = 1, 1, 7, 37, 266, 2431, ... . The solutions of Eqs. ()5.35|) and 
()5.36p are: q{x) = and v{x) = 0. It corresponds to 



e 2-''^ 



^ . ^(l-v/(2-at)-2A)a . 



in the boson formalism. 
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These examples show that any combinatorial structure which can be de- 
scribed by a Sheffer-type exponential generating function can be cast in bo- 
son language. This gives rise to a large number of formulas of the above type 
which put them in a quantum mechanical setting. 



Chapter 6 

MISCELLANY: APPLICATIONS 



Abstract 



We give three possible extensions of the above formahsm. First we show how 
this approach is modified when one considers deformations of the canonical boson 
algebra. Next we construct and analyze the generalized coherent states based on 
some number series which arise in the ordering problems for which we also furnish 
the solution to the moment problem. We end by using results of Chapter [S] to 
derive a substitution formula. 



6.1 Deformed bosons 

We solve the normal ordering problem for (A^A)"' where A and A^ are one 
mode deformed boson ladder operators ([A, A'*'] = [iV + 1] — [N]). The so- 
lution generalizes results known for canonical bosons (see Chapter E)). It 
involves combinatorial polynomials in the number operator N for which the 
generating function and explicit expressions are found. Simple deformations 
provide illustration of the method. 

Introduction 

We consider the general deformation of the boson algebra |Sol94j in the form 

[A, N] = A, 

[A\N] = -A\ (6.1) 
[A,A^ = [N+1]-[N]. 



Chapter 6. Miscellany: Applications 



60 



In the above A and are the (deformed) annihilation and creation op- 
erators, respectively, while the number operator counts particles. It is 
defined in the occupation number basis as N\n) = n\n) and commutes with 
A'^A (which is the consequence of the first two commutators in Eq. ()6.1|) ). 
Because of that in any representation of can be written in the form 

A'^A = [N], where [N] denotes an arbitrary function of A^, usually called 
the "box" function. Note that this is a generalization of the canonical boson 
algebra of Eqs. ()2.1|) and ()2.8|) described in Section (it is recovered for 
[N] = N + const). 

For general considerations we do not assume any realization of the number 
operator N and we treat it as an independent element of the algebra. More- 
over, we do not assume any particular form of the "box" function [A^]. Special 
cases, like the so(3) or so(2, 1) algebras, provide us with examples showing 
how such a general approach simplifies if an algebra and its realization are 
chosen. 

In the following we give the solution to the problem of normal ordering of a 
monomial (A^A)"' in deformed annihilation and creation operators. It is an 
extension of the problem for canonical pair [a, a^] = 1 described in Chapter 
El where we have considered Stirling numbers S{n, k) arising from 

n 

(a"^a)" = ^ S{n, h:){a)fa!' (for canonical bosons). (6.2) 

fc=i 

In the case of deformed bosons we cannot express the monomial {A^ A)"^ as 
a combination of normally ordered expressions in and A only. This was 
found for g-deformed bosons some time ago, |KK92j |Kat02j [KD95j |Sch03j . 
and recently for the R-deformed Heisenberg algebra related to the Calogero 
model, ^BN05] . The number operator A^ occurs in the final formulae because 
on commuting creation operators to the left we cannot get rid of A^ if it is 
assumed to be an independent element of the algebra. In general we can look 
for a solution of the form 

n 

{AAr = Y,Sn,k{N) {AfA\ (6.3) 

k=l 

where coefficients Sn,k{N) are functions of the number operator A^ and their 
shape depend on the box function [A^]. Following |KK92j we will call them 
operator-valued deformed Stirling numbers or just deformed Stirling numbers. 
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In the sequel we give a comprehensive analysis of this generalization. We 
shall give recurrences satisfied by Sn,k{N) of ()6.3|) and shall construct their 
generating functions. This will enable us to write down Sn,k{N) explicitly 
and to demonstrate how the method works on examples of simple Lie-type 
deformations of the canonical case. 

Recurrence relations 

One checks by induction that for each k > 1 the following relation holds 

[A^ A^] = {[N + k]- [N])A''~\ (6.4) 

Using this relation it is easy to check by induction that deformed Stirling 
numbers satisfy the recurrences^ 

Sn+l,k{N) = Sn,k-l{N) + {[N] - [iV - k])Sn,k{N) foT 1< k < u (6.5) 

with initial values 

5„,i(iV) = {[N] -[N- l]r-\ 5„,„(iV) = 1. (6.6) 

The proof can be deduced from the equalities 

ji+i 



k=l 



J2^n,kiN) (A^)' A'^At A'^5^5„,fc(iV) iA^f{A^A+[N + k]-[N])A' 



k=l k=l 

^ J2Sn,k-iiN) {A^)'A' + -[N- k])Sn,kiN) {A^)'A'. 

k=2 k=l 

It can be shown that each Sn,k{N) is a homogeneous polynomial of order 
n — k in variables [N], ...,[N — k]. For a polynomial box function one obtains 
<Sn,k{N) as a polynomial in A^. As we shall show, a simple example of the 
latter is the deformation given by the so{3) or so{2, 1) Lie algebras. 
Note that for the "box" function [N] = N + const {i.e., for the canonical 
algebra) we get the conventional Stirling numbers of Chapter 01 

^ The recurrence relation Eq. H6.5|l holds for all n and k if one puts the following "bound- 
ary conditions": Sij{N) = for « = or j = or i < j, except Sofl{N) = 1. 
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Generating functions and general expressions 

We define the set of ordinary generating functions of polynomials Sn,k{N), 
for A; > 1, in the form 

oo 

Pk{x,N) ■.= Y,Sn,k{N) x\ (6.7) 

n=k 

The initial conditions ()6.6|1 for k = 1 give 

Using the recurrencies ()6.5p completed with ()6.6|) one finds the relation 

whose proof is provided by the equalities 

Pfc(x, N)^f2 iSn-i,k-i{N) + {[N]-[N- k]) Sn-iAN)) 



n=k 

OD 



= xJ2^n-i,k-i{N) + {[N] - [AT - k]) xJ2Sn-iAN) 

n=k n=k 

= X Pk-i{x, N) + {[N]-[N- k]) X Pk{x, N). 

The expressions ()6.8|) and ()6.9|) give explicit formula for the ordinary gener- 
ating function 

For the canonical algebra it results in the ordinary generating function for 
Stirling numbers S{n, k): 

"^'^^^ = {i-x){i-2x).....{i-kxy 

Explicit knowledge of the ordinary generating functions ()6.1()|1 enables us to 
find the Sn,k{N) in a compact form. As a rational function of x it can be 
expressed as a sum of partial fractions 



Chapter 6. Miscellany: Applications 



63 



where 



Or = 7 ^. (6.13) 



j=l,j^r I [N]-[N-r] 



From the definition ()6.7|) we have that Sn,k{N) is the coefficient multiplying 
x"" in the formal Taylor expansion of Pk{x,N). Expanding the fractions in 
above equations and collecting the terms we get Sn,k{N) = Ylt=i^r{[N] — 
[N — r])"~'^. Finally, using (j6.13p . we arrive at 



\n-l 



where monotonicity of [A^] was assumed. 

For the canonical case it yields the conventional Stirling numbers S{n, k) = 
fE-=i see Eq.^. 

The first four families of general deformed Stirling numbers (polynomials) 
defined by ()6.3p read 



5i,i(iV) 


= 1, 






52,1 (iV) 


= [N]- 






52,2 (iV) 


= 1, 






53,1 (iV) 


= m - 






53,2 (iV) 


= m] - 


-[N-1]- 


[N 


53,3 (iV) 


= 1, 






54,1 (iV) 


= m - 


- [N - l])^ 




54,2 (iV) 


= 3[Ar]2 


-3[Ar][Ar- 


-1] 


54,3 (iV) 


= 3[iV] - 


-[N-l]- 


[N 


54,4(iV) 


= 1. 







2], 



3[A^] [N -2] + [N ~ 1]2 + [AT - 2] [A^ - 1] + [A^ - 2f 
2]-[iV-3], 



Examples of simple deformations 



If we fix the "box" function as [A^] = =F— 2 — then ()6.1|) become structural 
relations of so(3) and so{2, 1) algebras, respectively. The first four families 
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of polynomials, satisfying S^j. ' ( 
so(3) 

S,,i{N) = 1 

S2AN) = -N + l 
S2AN) = 1 

53,l(iV) = (iV-l)2 

SsAN) = -3iV + 4 

S3,3{N) = 1 

54,i(iV) = -(iV-l)3 
S4,31n) = -6N + 10 

sUn) = 1 



= (-ir-''5::f (iV),are 
so(2,l) 

Si,i{N) = 1 
52,i(iV) = iV-l 

52,2(iV) = 1 

Ss,i{N) = {N - ly 
53,2 (iV) = 3iV - 4 

53,3(iV) = 1 

54,1 (iV) = {N-lf 
54,2(A^) = 7iV2 _ i9Ar + 13 
S^^sIn) = 6A^ - 10 

sUn) = 1 



Up to now we have been considering A, A^ , and N as independent elements 
of the algebra. Choosing the representation and expressing in terms of 
A and A"^ we arrive at expressions like ()6.3|) which can be normally ordered 
further. We may consider a realization of the so{2, 1) in terms of canonical 
operators a and 

A = ^-=aa A^ = ^^a^a^ iV = - ( a'^a + ^ ) . (6.15) 

The solution of the normal ordering problem for ((a^)^a^)" was already con- 
sidered in Section |4. 7. H and leads to 

2n 

((at)2a2)" = ^52,2(n,fc)atV, (6.16) 

k=2 

where S2,2{n,k) are generalized Stirling numbers for which we have (see 
Eqs. (11121) and (jOEI)) 
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S2,2in,k) = ^E(-l)'Q[P(P-l)r (6-17) 
= E(-iy('^)5(2r.-/,fc). (6.18) 

Terms in with k even are directly expressible by A and A'^ . If k is odd 

then the factor a^a may be commuted to the left and we get (the convention 
of Ea. dOnil is used) 



((at)2a2)" = {S2,2{n, 2k) + S2,2{n, 2k + l){a^a - 2k)) 

k=l 

Finally, we get 

n 

(AU)" = J2 ('52,2(«, 2k) + S2,2{n, 2k + 1){2N -2k- 1/2)) A^'^A'^ 



k=l 



This indicates that the general solution involving higher order polynomials in 
can be simplified when a particular realization of the algebra is postulated. 
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6.2 Generalized coherent states 



We construct and analyze a family of coherent states built on sequences of 
integers originating from the solution of the boson normal ordering problem 
investigated in Section 14.7.11 These sequences generalize the conventional 
combinatorial Bell numbers and are shown to be moments of positive func- 
tions. Consequently, the resulting coherent states automatically satisfy the 
resolution of unity condition. In addition they display such non-classical 
fluctuation properties as super-Poissonian statistics and squeezing. 



Introduction 

Since their introduction in quantum optics many generalizations of standard 
coherent states have been proposed (see Appendix . The main purpose of 
such generalizations is to account for a full description of interacting quantum 
systems. The conventional coherent states provide a correct description of a 
typical non-interacting system, the harmonic oscillator. One formal approach 
to this problem is to redefine the standard boson creation and annihilation 
a operators, satisfying [a, a^] = 1, to A = af{a^a), where the function f{N), 
N = a^a, is chosen to adequately describe the interacting problem. Any 
deviation of f{x) from f{x) = const describes a non-linearity in the system. 
This amounts to introducing the modified (deformed) commutation relations 
Pol94j |MMHZ97j |MM98j (see also Section lO) 

[A, A^] = [N + 1] - [N], (6.19) 

where the "box" function [A^] is defined as [A^] = Np{N) > 0. Such a way of 
generalizing the boson commutator naturally leads to generalized "nonlinear" 
coherent states in the form {[n]\ = [0][1] . . .[n], [0] = 1) 




''n\l 



n=0 

which are eigenstates of the "deformed" boson annihilation operator A 

A\z)=z\z). (6.21 



It is worth pointing out that such nonlinear coherent states have been suc- 
cessfully applied to a large class of physical problems in quantum optics 
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|MV96j [kVD Ol]. The comprehensive treatment of coherent states of the 
form of Eq. (IFT^ can be found in [Mi] [MM^tZQT] |MM98j . 
An essential ingredient in the definition of coherent states is the completeness 
property (or the resolution of unity condition) KS85J [Kla63aj [Kla63bJ [KPSOlj. 
A guideline for the construction of coherent states in general has been put 
forward in |Kla63aj as a minimal set of conditions. Apart from the conditions 
of normalizability and continuity in the complex label this set reduces to 
satisfaction of the resolution of unity condition. This implies the existence 
of a positive function VF(|2;p) satisfying [KPSOlj 



d^z \z)W{\z\^){z\ 



^\n){n\ 



(6.22) 



n=0 



which reflects the completeness of the set {\z)}. 

In Eq. ()6.22|) / is the unit operator and \n) is a complete set of orthonormal 
eigenvectors. In a general approach one chooses strictly positive parameters 
p{n), n = 0,1, . . . such that the state \z) which is normalized, {z\z) = 1, is 
given by 



\z) = Af-'/\\z\^)Y: 



n=0 



(6.23) 



with normalization 



n=0 



(6.24) 



which we assume here to be a convergent series in for all z G C. In view of 
Eqs. fl6.19|) and ()6.20|) this corresponds to p{n) = [n]\ or [n] = p{n)/ p(n — 1) 
for n = 1,2, . . .. 

Condition ()6.22|) can be shown to be equivalent to the following infinite set 
of equations {KPSOlj : 



X 



W{x) 
Af{x) 



dx = p{n). 



n 



0,1, 



(6.25) 



which is a Stieltjes moment problem for Vr(x) = tt . 

Recently considerable progress was made in finding explicit solutions of Eq. ()6.25jl 
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for a large set of p{n)^s, generalizing the conventional choice \z)c for which 
Pc{n) = n\ with Mc{x) = (see Refs. |KPSnij [P^ |SPKnij ^MU] 
and references therein), thereby extending the known families of coherent 
states. This progress was facilitated by the observation that when the mo- 
ments form certain combinatorial sequences a solution of the associated 
Stieltjes moment problem may be obtained explicitly |PS01j |BPS03aj [BPS04j . 

In the following we make contact with the combinatorial sequences ap- 
pearing in the solution of the boson normal ordering problem considered in 
Section l4.7.1l These sequences have the very desirable property of being mo- 
ments of Stieltjes- type measures and so automatically fulfill the resolution of 
unity requirement which is a consequence of the Dobihski-type relations. It is 
therefore natural to use these sequences for the coherent states construction, 
thereby providing a link between the quantum states and the combinatorial 
structures. 

Moment problem 

We recall the definition of the generalized Stirling numbers Sr,i{n, k) of Sec- 
tion HTiU defined through (n, r > integers): 

n 

[(at)"a]" = (at)"("-^)^S,,l(r^,fc)(a^)'=a^ (6.26) 
fe=i 

as well as generalized Bell numbers Br^i{n) 

n 

Br,i{n) =Y,Sr,i{n,k). (6.27) 



k=l 



For both Sr,i{n,k) and Br^i{n) exact and explicit formulas have been given 
in Section 14.7.11 

For series Br^i{n) a convenient infinite series representation may be given by 
the Dobihski-type relations (see Eas. lHTni) . and (nmUl ) 

1 r(n + ^) 

ft,,(„)^L_L5^_^^, r >l, (6,28) 

fc=0 ^ r-l' 



and 



k\ 

k=0 
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From Eqs. (j6.26j) and (j6.'27|) one directly sees that Br^i{n) are integers. Ex- 
plicitly: = 1,2,5,15,52,203,...; B2,i{n) = l', 3, 13, 73, 501, 4051 . . .; 
53,i(n) = 1^ 4, 25, 211, 2236, 28471, . . .; 54,i(n) = 1, 5, 41, 465, 6721, 117941 . . . 
etc. 

It is essential for our purposes to observe that the integer Br^iin + 1), n = 
0, 1, ... is the n-th moment of a positive function Wr^i{x) on the positive half- 
axis. For r = 1 we see from Eq. ()6.29|) that Bi i{n + 1) is the n-th moment of 
a discrete distribution Wi^i{x) located at positive integers, a so-called Dirac 
comb: 

POO 

+ 1) = X 
Jo 

For every r > 1 a continuous distribution Wr^i{x) will be obtained by excising 
(r — 1)"" r(n + ^3y) from Eq. ()6.28|) . performing the inverse Mellin transform 
on it and inserting the result back in the sum of Eq. ()6.28|) . see pCPSOlJ [SixOlj 
for details. (Note that -Br-,i(0) = r = 2, 3, ... is no longer integral). In 
this way we obtain 

/■oo 

Br,iin + 1)= x'^Wr,iix) dx, (6.31) 
Jo 

which yields for r = 2, 3, 4: 

W2,i{x) = e-"-^v^ /i(2v^), (6.32) 

W„(.) = -^-e-.-' (-^ „F,(-.-- -) + ^ „F.(-, 2; -) j (6.33) 

4 2 ^ ,2 4 5 X , ^,2, ^ ,4 5 „ X A 
3'-' „^'3(3.3.3;5^)+>rr(-)x„F,(-.-.2;-)j. 

In Eqs. (j6.32|) . (j6.33p and (j6.34|) Iu{y) and o-^p(- • • ; 2/) are modified Bessel and 
hypergeometric functions, respectively. Other Wr^i{x) > for r > 4 can be 
generated by essentially the same procedure. 

In Figure IH?T] we display the weight functions Wr^i{x) for r = 1 . . . 4; all of 
them are normalized to one. In the inset the height of the vertical line at 
X = k symbolizes the strength of the delta function 6{x — k), see Eq. (l6.30|l . 
For further properties of Wi^i{x) and more generally of Wr^r{x) associated 
with Eq. fjFTTfj) . see |HPS04j . ' 



^ oo 



(k 

k=i ^ 



6{x — k) 
1)! 



dx. 



(6.30) 
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0.25 




X 



Fig. 6.1: The weight functions Wr^i{x), (x = in the resolution of unity 
for r = 2,3,4 (continuous curves) and for r = 1, a Dirac's comb 
(in the inset), function of x. 



Construction and properties 

A comparison of Eqs. (j6.23p . (j6.24|) and (j6.31|) indicates that the normaUzed 
states defined through p{n) = Br^i{n + 1) as 

oo „ 

\z), = K^'Wz\^) y , |n), (6.35) 



with normahzation 

oo 

X 

n=0 



automatically satisfy the resolution of unity condition of Eq. ()6.22|) . since for 

„ oo 
/ (fz \z)rWrl{\z\^)r{z\ = ^ = (6-37) 
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Fig. 6.2: Mandel parameters Qr{x), forr = 1 ... 4, as a function ofx= \z 
see Eq. X^rm . 



Note that Eq. ()6.35|) is equivalent to Eq. ()6.20|) with the definition [n\r = 
5,,i(n + l)/5,,i(n),n = 0,1,2,.... 

Having satisfied the completeness condition with the functions Wr^i{x), r = 
1,2,... we now proceed to examine the quantum-optical fiuctuation proper- 
ties of the states \z)r. From now on we consider the |n)'s to be eigenfunctions 
of the boson number operator = a^a, i.e. N\n) = n\n). The Mandel pa- 
rameter jKPSOlj 



allows one to distinguish between the sub-Poissonian (antibunching effect, 
Qr < 0) and super-Poissonian (bunching effect, Qr > 0) statistics of the 
beam. In Figure we display Qr{x) for r = 1 . . .4. It can be seen that all 
the states \z)r in question are super-Poissonian in nature, with the deviation 
from Qr = 0, which characterizes the conventional coherent states, diminish- 
ing for r increasing. 




(6.38) 
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1.8 




0.6 - 







2 



3 



4 



Re(z) 



Fig. 6.3: The squeezing parameters of Eqs. \6.39}) and il6'.40)) for tlie coor- 
dinate Q (three upper curves) and for the momentum P (three 
lower curves) for different r, as a function of Re{z), for r = 1,2, 3. 



In Figure we show the behavior of 



which are the measures of squeezing in the coordinate and momentum quadra- 
tures respectively. In the display we have chosen the section along Re{z). All 
the states l^;),. are squeezed in the momentum P and dilated in the coordinate 
Q. The degree of squeezing and dilation diminishes with increasing r. By 
introducing the imaginary part in z the curves of Sq{z) and Sp{z) smoothly 
transform into one another, with the identification Sqlia) = Sp{a) and 
Sp{ia) = Sgi^a) for any positive a. 

In Figure we show the signal-to-quantum noise ratio |Yue76j relative to 




(6.39) 



and 




(6.40) 
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4[c(2;|A^|z)c] = 4|2;p, its value in conventional coherent states; i.e the quantity 
o-j. = (Tj. — 4:[c{z\N\z)c], where 



'{z\Q\z)r 

(AQ)2 



(6.41) 



with {AQY = r{z\Q'^\r)r — {r{z\Q\r)r)'^. Again only the section Re{z) is 
shown. We conclude from Figure IHTil that the states \z)r are more "noisy" 
than the standard coherent states with pdn) = n\ . 
In Figure we give the metric factors 



Afr'i 



X 



Mr 



(6.42) 



which describe the geometrical properties of embedding the surface of coher- 
ent states in Hilbert space, or equivalently a measure of a distortion of the 
complex plane induced by the coherent states [KPSOlj . Here, as far as r is 
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concerned, the state \z)i appears to be most distant from the \z)c coherent 
states for which Uc = 1- 

Remarks 

The use of sequences Br^i{n) to construct coherent states is not hmited to the 
case exemphfied by Eq. ()6.35|) . In fact, any sequence of the form Br^i{n 
p = 0, 1, ... will also define a set of coherent states, as then their respective 
weight functions will be V^{x) = > 0. Will the physical prop- 

erties of coherent states defined with Pp{n) = Br^iin + p) depend sensitively 
on p ? A case in point is that of p = for which qualitative differences from 
Figure 16.21 (p = 1) appear. In Figure 16.61 we present the Mandel parame- 
ter for these states. Whereas for r = 1 the state is still super-Poissonian, 
for r = 2,3,4 one observes novel behavior, namely a crossover from sub- to 
super-Poissonian statistics for finite values of x. Also, as indicated in Fig- 
ure 16.71 we note a cross-over between squeezing and dilating behavior for 
different r. These curious features merit further investigation. 
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Fig. 6.7: Squeezing parameters for the states with p(n) = Bj.^i{n) for r = 
1, 2, 3, as a function of x = Iz]"^. The subscripts P and Q refer to 
momentum and coordinate variables respectively. 
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6.3 Substitution theorem 

In Chapter El we have demonstrated how to find the normally ordered form 
of the exponential of the operator linear in the annihilation or creation op- 
erator. This representation is of fundamental meaning in the coherent state 
representation used in physics. 

As a demonstration of possible other utility of such formulas we show the di- 
rect proof of the substitution theorem. More specifically we will demonstrate 
that for any function F{x) the following equality holds 

^x(,ix)D+v(x))p^^^ = ^(A, x) ■ F(T(A, x)) (6.43) 

where functions T(A, x) and G{X,x) may be found from the following equa- 
tions 

^l2^ = q(T{X,x)) , T{0,x)=x, (6.44) 



dX 



viT{X,x))-GiX,x) , G{0,x) = l. (6.45) 



First observe from Eq. ()6.43j] that the action of e''^'-''*^^'''^"''^*^^^'' on a function 
F{x) amounts to: 

a) change of argument x T(A, x) in F{x) which is in fact a substitution; 

b) multiplication by a prefactor G{X, x) which we call a prefunction. 

We also see from Eq. ()6.45p that G{X,x) = 1 for v{x) = 0. Finally, note 
that e'*''-'^^"^)^"'"^*^"^^) with A real generates an abelian, one-parameter group, 
implemented by Eq. (j6.43j) : this gives the following group composition law for 
T(A, x) and G(A, x): 

TiX + e,x) = mT(A,x)), 

G{X + 9,x) = G{X,x) ■G{e,T{X,x)). ^ ^ ^ 

Now we proceed to the proof of Eqs. ()6.43|) - ()6.45|) . Using the multiplication 
X and derivative D representation we may rewrite Eq. ()5.28|) as 



with accompanying differential Eqs. (j5.25|) and (j5.26p for functions T(A,x) 
and G{X,x) which are exactly the same as Eqs. (j6.44|) and (|6.45p . 



Chapter 6. Miscellany: Applications 



77 



Now, recalling the Taylor formula e"^F{x) = F{x + a) we may generalize it 
to 

: e"(^)^ : = + (6.48) 

It can be seen from the power series expansion of the exponential. 
This allows to write 

= G(A,X) : e[^(^'^)-^l^ : 

G{X,X) F{x+[T{X,x)-x]) 
= G{X,x) ■ F{T{X,x)). 

This completes the proof. 

Here we list several examples illustrating Eqs. ()6.43p - ()6.45|) for some choices of 
q{x) and v{x). Since Eqs. (l6.44p and (|6.45|) are first order linear differential 
equations we shall simply write down their solutions without dwelling on 
details. First we treat the case of v{x) = 0, which implies g{X, x) = 1: 

Ex.1: 

q{x) = X, T(A,x) = xe'^ 

which gives exp (Ax^) F{x) = F(xe^), a well known illustration of the 
Euler dilation operator exp (Ax^). 

Ex.2: 

X 



g(x) = x*", r > 1, T(A,x) 



(1 - A(r - l)x'-^)^-' 



The above examples were considered in the literature |D()TV97j [Lannnj |BPSn8bj |BPSn3cj . 
We shall go on to examples with v{x) ^ leading to nontrivial prefunctions: 



Ex.3: 



g(x) = 1 and v[x) — arbitrary, 
T(A, x) = X + A, 

A 

/ vix^u) 

g{X,x) = eo 
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Ex.4: 



q{x) = X, and v{x) = x^ 
T(A, x) = xe^, 



g{\,x) = exp 



y (e" - 1) 



Ex.5: 



q{x) = r > 1 and v{x) = x'^ 
T{X,x) 



g{X,x) = exp 



(1 - A(r-l)a;"-i)'-^ 

-r+l / 2 



X 



1 — r 



(l-A(r-l)x"-^) 



Closer inspection of the above examples (or at any other example which 
the reader may easily construct) indicates that from the analytical point of 
view the substitution theorem should be supplemented by some additional 
assumptions, like restrictions on A. In general we can say that it can be 
valid only locally |DPS"'"n4j . However, we note that a formulation in the 
language of formal power series does not require such detailed analysis and 
the resulting restrictions may be checked afterwords. 

We finally note that the substitution theorem can be used to find the normally 
ordered form of an exponential operator linear in a or a"!" (as in Chapter E}. 
This idea has been exploited in |BHP+05] . 



Chapter 7 

CONCLUSIONS 



In this work we have considered the boson normal ordering problem for pow- 
ers and exponentials of two wide classes of operators. The first one consists of 
boson strings and more generally homogenous polynomials, while the second 
one treats operators linear in one of the creation or annihilation operators. 
We have used the methods of advanced combinatorial analysis to obtain a 
thorough understanding and efficient use of the proposed formalism. In all 
cases we provided closed form expressions, generating functions, recurrences 
etc. The analysis was based on the Dobihski-type relations and the umbral 
calculus methods. In general, the combinatorial analysis is shown to be an 
effective and flexible tool in this kind of problem. 

We also provided a wealth of examples and pointed out possible applications. 
The advantages may be neatly seen from the coherent state perspective {e.g. 
we may use it for construction of the phase space pictures of quantum me- 
chanics). We may also obtain solutions of the moment problem, enabling 
us to construct new families of generalized coherent states. Moreover appli- 
cation to the operator calculus is noted and exemplified in the substitution 
theorem. We also observe that the normal ordering problems for deformed 
algebras may be handled within that setting. 

These few remarks on the possible applications and extensions of the meth- 
ods used in this work (others may be found in the supplied references) yield 
a potentially wide field of possible future research. 

We would like to point out some other interesting features. We believe that 
similar combinatorial methods may be used to find the normally ordered 
forms of other classes of operators, such as the exponential of a general boson 
polynomial. Moreover, the extension to the multi-mode boson and fermion 
case should be susceptible to explicit analysis. This program would open a 
wide arena for immediate application to variety of physical models. 
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In conclusion, we would like to recall that the use of advanced combinatorics 
provides a wealth of interpretative tools for the problems under consider- 
ation. Combinatorial objects may be interpreted in terms of graphs, rook 
polynomials, partitions, correlations etc. We have not touched upon this 
aspect in this work but we consider it as a very promising one. 



Appendix A 

COHERENT STATES 



Here we define and review some properties of cofierent states |KS85j |ZFR90j . 
We define coherent states l^;) as the eigenstates of the annihilation operator 

a\z) = z\z), (A.l) 

where 2; G C is a complex number. They can be written explicitly as 



00 

1 l,|2 



z) = e-^\^\' ^^|n). (A.2) 



Ov ^ 



These states are normalized but not orthogonal. The coherent states over- 
lapping factor is 

{z\z') = e^*^'-f l-V-lkP^ (A.3) 

The set of coherent states {1^) : z E C} constitute an overcomplete basis in 
the Hilbert space Ti. The following resolution of unity property holds 

- / \z){z\ d^z = l. (A.4) 

Eq. ()A.2jl may be rewritten in the form (see Eqs. ()2.(')|l ) 

1^) = e-^e""'|0), (A.5) 
or with the use of Eq. ()2.13|1 as 

\z) = e'''^-'*^\0). (A.6) 

The operator D called a displacement operator and therefore 

coherent states are sometimes called displaced vacuum. This group property 
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may be taken as the definition of the coherent states also for other (than 

Heisenberg-Weyl) groups, see |Per86j . 

Let us define two self-adjoint operators Q and P by 



Q 



f ' (A.7) 

^ a' — a 
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They satisfy the Heisenberg-Weyl commutation relation [Q, P] = i. These 
operators may be interpreted as the position and momentum operators for 
the quantum particle in the harmonic oscillator potential, in quantum optics 
they play the role of the field quadratures and are also used in the phase space 
formulation of Quantum Mechanics. It can be shown that the coherent states 
are the minimum uncertainty states for operators Q and P, i. e. 

A|,)g-A|,)P = i, (A.8) 

where A^A = a/ {ip\{A — {ip\A\ip)Y\ip) is the uncertainty of the operator A. 
Moreover they are the only states if one additionally imposes the condition 
A.\z)Q = ^\z)P (otherwise the family of squeezed states is obtained) which 
serves as another possible definition of coherent states. 
Finally we mention that the resolution of unity of Eq. ()A.4|) together with 
the continuity of the mapping z ^ \z) are sometimes taken as the mini- 
mum requirements for the coherent states. However, this definition is not 
unique and leads to other families of so called generalized coherent states 
|K1afi3aj [Kk63b] [K^ (see also Section ES)). 

We note that, due to their special features, coherent states are widely used 
in quantum optics [Gla 63^ KS(2I ^ ^^^^ other areas of physics [KS85j . 
In this text we especially exploit the property of Eq. ()A.1|) . It is because for 
an operator P(a, a^) which is in the normal form, P(a, a^) = Af [P(a, a^)] = 
: P(a,a^) :, its coherent state matrix elements may be readily written as 

{z\F{a,a^)\z') = {z\z') F{z',z*). (A.9) 

Also for the double dot operation it immediately yields 

{z\ : G{a, a^) : \z') = {z\z') G{z', z*). (A.IO) 
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Unfortunately for the general operator none of these formulae hold. Never- 
theless, there is a very useful property which is true; that is, if for an arbitrary 
operator F{a, a^) we have 

{z\F{a,a^)\z') = {z\z') G{z*,z') (A.ll) 

then the normally ordered form of F{a, a"!") is given by 

A/'[F(a,a^)] = : G{a\a) : . (A.12) 

For other properties of coherent states we refer to |KS85j|ZFR90j . 



Appendix B 



FORMAL POWER SERIES. UMBRAL 
CALCULUS 



Here we recall the basic definitions and theorems concerning the formal power 
series calculus. As an illustration some topics of the umbral calculus are re- 
viewed. For a detailed discussion see |JNiv69j . |Com74j . jRio84j . |Wil94j . 
|C4KF94j . [EHOSl, jHomHU, |H,ai65j . \HR7H\ . \Bnc9H\ . 



Formal power series 

Suppose we are given a series of numbers (/n)^o- define a formal power 
series in indeterminate X clS 

n=0 

The set of formal power series constitutes a ring when the following opera- 
tions are imposed 

• addition : 

oo ^ oo ^ oo 

F{x) + G{x) = E /"li + E ^-li = + 3n)'^_ (B.2) 

71=0 ?i=0 n=0 

• multiplication (Cauchy product rule): 

oo jjOO jjOon^\ 

n=0 ■ n=0 ■ n=0 fc=0 ^ ^ 
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Both operations have an inverses. The inverse of the series F{x) with respect 
to addition is the series with negative coefficients —F{x) = ~f"'^- The 

(unique) multiphcative inverse on the other hand is well defined only when 
/o 7^ and may be given recursively, see e.g. |,Wil94ij . 
One can also define the substitution of formal power series by 

F{G{x)) = f^^^ (B.4) 

n=0 

when = in the expansion G{x) = J2'^=o dn'n' ■ Explicit expression for the 
coefficients is given by the Faa di Bruno formula |Com74j |Ald01j . A series 
F{x) may be shown to have a compositional inverse F~^{x) iff /o = and 

Other definitions and properties can be further given. Here is the example 
of the derivative and integral of the formal power series (note that these 
operators act like shift operators on the sequence (/n)^o) 



DF{x) = F\x) = Y: fnj^, = E /..i^^ (B.5) 



fn+l 

n=l ' ?i=0 



Fix) dx = Y^ fnr—lv = E /-I- (B.6) 
^ (n + 1)! ^ n\ 

Observe that these definitions mirror the corresponding operations on ana- 
lytic functions. Nevertheless the limiting operations are not needed in these 
definitions. The ring of formal power series can be also be given the structure 
of a complete metric space (see e.g. |FS05j . |Rom84j ). 

We note that {fn)^=o may be also a sequence of polynomials, functions or 
anything else. 

The use of power series is well suited in the context of generating functions. 
The series as in Eq. ()B.l|) is called the exponential generating function (be- 
cause of the factor l/n\ in the expansion). Functions G{x) = Yl^=o9nX^ are 
called ordinary generating functions. The use of generating functions is es- 
pecially useful in solving recurrences and enumerating combinatorial objects, 
see eg. | Wil94j [HKPQij jF^ . 

All these definitions and properties may be naturally extended to the multi- 
variable case. 

For a systematic description of the formal power series and applications see 
|Nw69l . \(]om1i\ . jEioHll, |Wil94| . |(jKP94| . (EHOSI, |l^.om84| . 
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Umbral calculus. ShefFer-type polynomials 

Subject of the umbral calculus is the study of a Sheffer A-type zero poly- 
nomials, called here briefly Sheffer- type |She39j . |Rai65j . Without going 
into details we recall some basic theorems and definitions |R,om84j . |R,R,78j . 

Suppose we have a polynomial sequence s„(x), n = 0,1,2,... being 
the polynomial of degree n). It is called of a Sheffer-type if it possesses an 
exponential generating function of the form 

GiX,x) = = A(A)e^'^(^). (B.7) 

n=0 

for some (possibly formal) functions A{X) and B{X) such that -B(O) = 0, 
-B'(O) 7^ and A{0) ^ 0. When B[X) = 1 it is called .4^0^/6 sequence for 
A[X). For A(A) = 1 it is known as associated sequence for B{X). 
Yet another definition of the Sheffer-type sequences can be given through 
their lowering and rising operators, i.e. the polynomial sequence s„(x) is of 
Sheffer-type iff there exist some functions f{x) and g{x) (possibly formal) 
such that /(O) = 0, /'(O) ^ and ^(0) ^ which satisfy 

fiD)s^{x) = ns„„i(x), (B.8) 



^ (x) = Sn+i{x). (B.9) 



These two definitions describe the Sheffer-type sequence uniquely and the 
correspondence is given by 

A(x)^r\x), (B.IO) 

Many curious properties of the Sheffer-type polynomials can be worked out. 
We quote only one of them, called the Sheffer identity., to show the simplicity 
of formal manipulations. Using definition Eq. ()B.7|l we have 



ys^{x + y)^ = A(A)e(^+^)^W = A(A)e^^We^^W 

n=0 

n=0 ■ n=0 
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By the Cauchy product rule Eq. ()B.3|) we obtain the Sheffer identity: 

Sn{x + y) = Y1 (^Pkiy)Sn-k{x) (B.12) 

where Pn{x) is the associated sequence for -B(A). Observe that for associated 
sequences {A{X) = 1) this property generahzes the binomial identity and 
therefore sometimes they are called of binomial type. 

The Sheffer-type polynomials constitute a basis in the space of polynomials 
and through that property any formal power series can be developed in that 
basis also. On the other hand, formal power series can be treated as func- 
tional on the space of polynomials. Investigation of this connection is also 
the subject of umbral calculus in which the Sheffer-type polynomials play a 
prominent role |Rom84j . 



We emphasize that all the functions and series here can be formal and then 
the operations are meant in the sense of the previous Section. 
Thorough discussion of the umbral calculus, Sheffer-type polynomials and 
applications can be found in |Rom84j |RR78j |Buc98j |Rai65j . 
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